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I 6.  -^lb»»rOCt 

^■To  assist  in  the  evaluation  of  proposed  structural  inspection  programs  for 
commercial  jet  transport  aircraft,  a logic  was  developed  to  simulate  structural 
defects,  failures,  and  inspections.  This  logic  was  incorporated  in  a computer 
program  entitled  Structural  Area  Inspection  Frequency  Evaluation  (SAIFE).  With  the 
objective  of  quantifying  the  evaluation  process  currently  used  to  establish  and 
modify  inspection  intervals,  SAIFE  accounts  for  the  following  factors:  (1)  aircraft 
design  analysis;  (2)  fatigue  testing;  (3)  production,  service,  and  corrosion  de- 
fects; (4)  probability  of  crack  or  corrosion  detection;  and  (5)  aircraft  modifica- 
tion economics.  As  a five-volume  document,  this  report  covers  the  initial  contract 
effort  plus  a subsequent  parametric  analysis  as  follows:  Volume^  (antitiad  Execg- 
tAve^-Summary}-  presents  the  SAIFE  logic  and  documents  the  methodology  for  the  de- 
cision-making processes  in  the  simulation  logic.  Volume,  Ij  (-enfrit  1-ad  lies c r.i-p4men--»f 
>Sinwlatlan  Irogmtr)' details  the  SAIFE  simulation  logic,  presents  the  background  data 
for  the  analytical  functions  and  decision-making  processes,  and  includes  data  for  a 
typical  simulation  problem.  Volume  |I I tl cd-Etewcms Lra troir -Inptrty- inspection? 

presents  data  tabulations  derived  from  historical  trends  and 
design  input  data  for  a SAIFE  demonstration  problem.  As  the  user's  manual  for  the 
SAIFE  computer  program,  VolumeUVycontains  detailed  computer  logic  flow  diagrams  and 
a complete  listing  of  the  program  which  is  written  in  SIMSCRIPT  1 1.5.  Volume  V 
(■an  titled-  Resutt* of  -HodcT  Demon's  trartlen)''  presents  the  results  of  the  program  appli- 
cation to  a hypothetical  aircraft  and  compares  these  results  with  the  service  ex- 
perience of  operational  aircraft. 

17,  K *y  Words  ^ 

Aircraft  experience  simulation;  aircr'aftl 


J 


it 


inspection  intervals;  aircraft  hazard 
rate  and  reliability;  production,  servici 
and  corrosion  defects;  crack  or  corrosio! 


18.  Dlttr.bution  5tat.nt.nt 


' 


Document  is  available  to  the  public 
through  the  National  Technical 
Information  Service,  Springfield, 
Virginia  22151 . 


19.  Security  Clou  /,  {of  ihii  report) 

20.  Security  Class!/,  (of  this  pog*) 

21*  No.  of  Pog«i  1 

22.  Pric* 

UNCLASSIFIED 

UNCLASSIFIED 

280 

/ 

Form  DOT  F 1700.7  (8-72) 


Reproduction  of  completed  poge  authorised 

?//  -/ « 


w 


/7, 


METRIC  CONVERSION  FACTORS 


TABLE  OF  CONTENTS 


INTRODUCTION  1 

SYSTEM  DESCRIPTION  1 

GENERAL  PROGRAM  DESCRIPTION  3 

PROGRAM  INPUT  3 

AIRCRAFT  DATA 4 

LONG  LIST  DATA 7 

ELEMENT  DATA 7 

FORMAT  SPECIFICATIONS 11 

SAMPLE  INPUT  1 11 

PROGRAM  OUTPUT 14 

STANDARD  OUTPUT  ELEMENT  DATA 14 

STANDARD  OUTPUT  SUMMARY  DATA 28 

LONG  LIST  DATA 28 


APPENDIX  A - DETAILED  PROGRAM  DESCRIPTION 34 

APPENDIX  B - MODIFIED  PROGRAM  SOURCE  LISTING.  ...  92 


ACf'-'IS'?1:?!  tnr 

fills  v. Section 

0,10  B.:(i  Sno’lofi 

I'KAiWiNCtO  □ 

JlSU.'iOI'.'N  ■ 

I BY 

oisinB'jiffii/WJi'fVjiiiii 


INTRODUCTION 


The  Structural  Area  Inspection  Frequency  Evaluation  (SAIFE)  computer 
simulation  was  developed  for  the  FAA  under  contract  D0T-FA74WA-3493  by 
Technology  Inc.,  Dayton,  Ohio.  The  computer  program  simulates  the 
structural  inspection  program  of  an  airline  fleet.  SAIFE  accounts  for 
the  following  factors: 

1.  Aircraft  design  analysis 

2.  Component  and  full  scale  fatigue  testing 

3.  Production,  service  and  corrosion  defects 

4.  Probability  of  crack  or  corrosion  detection 

5.  Aircraft  component  repair  or  modification 

The  original  computer  model  has  been  extensively  modified  by  the  Engineering 
and  Manufacturing  Branch,  Flight  Standards  National  Field  Office. 
Modifications  have  occurred  in  program  input,  simulation  logic,  and 
program  output.  Major  modifications  were: 

1.  Replacement  of  the  two  crack  growth  rate  model  by  a four  crack 
growth  rate  model . 

2.  Introduction  of  sampling  at  the  D-level  inspection  in  that 
certain  aircraft  receive  an  Internal  inspection  based  upon  program  input 
and  program  logic. 

3.  Additional  logic  to  calculate  the  probability  of  aircraft 
failure  based  upon  crack  distribution. 

4.  Changes  to  inspection  reliability  curves  which  determine 
probability  of  defect  detection  at  a given  inspection  level. 

5.  Changes  to  logic  for  reducing  Inspection  intervals  and  scheduling 
special  inspections  based  on  crack  detection. 

During  the  program  modification  period  Technology  Inc.  provided  support 
under  an  additional  FAA  contract. 


SYSTEM  DESCRIPTION 

The  eight  blocks  in  Figure  1 represent  the  major  aspects  of  the  SAIFE 
simulation  logic.  Clock  1 accepts  input  data  for  the  aircraft  fleet  and 
for  each  structural  element  in  the  aircraft.  After  determining  whether 
element  modifications  are  required  because  of  the  fatigue  test  results 
in  Block  2,  Block  1 assigns  a fatigue  life  to  each  element  in  each 
aircraft.  Block  3 determines  whether  production,  service,  or  corrosion 
defects  will  occur;  if  it  is  determined  that  such  defects  will  occur. 
Block  3 predicts  the  times  when  they  will  occur.  After  comparing  the 
flight  loads  with  the  strength  of  each  element.  Block  4 predicts  the 
time  to  failure  for  each  element.  Block  5 conducts  the  periodic  inspec- 
tions of  each  element.  If  defects  are  detected,  Block  6 repairs  the 
element  and  assigns  it  a new  fatigue  life.  However,  if  an  existing 
defect  is  allowed  to  grow  until  element  failure.  Block  5 deletes  the 
aircraft  from  the  fleet.  Depending  on  the  magnitude  of  the  detected 
defects,  special  inspections  and  increased  inspection  frequencies  may  be 
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Figure  1.  Major  Aspects  of  SAIFE  Logic 
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called  for  in  Block  7 and  modifications  may  be  instituted  In  Block  8. 
When  all  the  aircraft  have  been  deleted  from  the  fleet  or  retired  from 
service,  the  simulation  is  complete. 


GENERAL  PROGRAM  DESCRIPTION 

The  SAIFE  program  is  a large,  complex  math  model  designed  to  simulate 
the  structural  performance  of  an  aircraft  fleet  and  the  effectiveness  of 
the  inspection  program  for  the  ait-craft  fleet.  The  aircraft  model  in 
the  computer  simulation  is  divided  into  structurally  significant  elements 
and  the  Inspection  program  for  each  element  is  defined.  Structural 
defects  are  classified  as  follows:  fatigue  and  corrosion  which  are 
wear-out  and  aging  phenomena;  production  or  design  defects;  and,  opera- 
tional or  maintenance  damage.  These  defects  and  the  inspection  program 
are  treated  as  probabilistic  phenomena  interacting  o */er  time.  SIMSCRIPT 
II. 5 (trademark.  Consolidated  Analysis  Centers,  Inc.,  Los  Angeles, 
California),  a computer  language  designed  for  discrete-event  simulation 
applications,  was  chosen  for  the  SAIFE  simulation. 

SAIFE  uses  ten  random  number  streams  to  generate  the  probabilistic 
phenomena  in  the  simulation.  Because  random  numbers  generation  is 
computer  word  length  dependent,  different  computer  systems  will  generate 
different  random  numbers  and  thus  different  output  for  a given  set  of 
input  data.  This  has  been  noticed  since  the  computer  program  has  been 
run  on  both  IBM  and  CDC  computers.  Results  may  be  duplicated  on  a given 
computer  since  the  compiler  initializes  the  seeds  for  the  random  number 
streams  to  the  same  value  before  each  run.  An  option  has  been  incor- 
porated in  the  computer  program  so  that  the  user  can  input  seeds  to 
evaluate  independent  trials  for  a given  set  of  input. 

Because  of  the  extensive  detail,  the  program  events  and  routines  are 
described  in  Appendix  A. 


PROGRAM  INPUT 

The  program  input  consists  of  three  parts.  The  first  part  consists  of 
input  variables  which  pertain  to  the  aircraft  type  under  consideration. 
These  variables  are  input  only  once  per  simulation  run  and  are  constant 
from  element  to  element.  If  the  user  desires  to  input  random  number 
seeds,  the  ten  seeds  are  input  after  the  aircraft  input.  The  second 
part  of  the  input  consists  of  long  list  data  if  standard  output  is  not 
used.  The  third  part  of  the  input  consists  of  input  variables  whose 
values  are  unique  to  each  element.  These  variables  must  be  input  in 
their  entirety  for  each  element  being  simulated. 
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AIRCRAFT  DATA.  The  input  variables  which  pertain  to  the  aircraft  type 
are  listed  and  described  below  in  the  order  in  which  they  are  read  in 
by  SAIFE. 

MODEL  (*}  - This  one-dimensional  alpha  array  of  size  two  Identifies 
the  aircraft  type  under  consideration.  The  total  length  of  this 
Identification  cannot  exceed  eight  characters. 

SIZE. OF. FLEET  - This  integer  variable  is  the  number  of  aircraft  in 
the  fleet  being  simulated.  The  output  format  requires  that  this 
variable  does  not  exceed  99999. 

USAGE. LIFE  - This  real  variable  is  the  service  life  in  flight  hours 
of  the  aircraft  being  simulated.  All  aircraft  in  the  fleet  must 
have  the  same  service  life.  The  output  format  requires  that  this 
variable  does  not  exceed  9999999. 

BEGIN. PRODUCTION  - This  real  variable  is  the  time  in  flight  hours 
rela t i ve  to  the  start  of  the  simulation  when  the  first  aircraft 
enters  service.  This  variable  in  conjunction  with  the  input  variable 
START. TEST  enables  the  user  to  start  the  fatigue  test  of  the  element 
before,  after,  or  at  the  same  time  the  first  aircraft  enters  service. 

PRODUCTION. TIME  - This  real  variable  defines  the  initial  aircraft 
production  rate.  It  is  the  time  In  flight  hours  between  aircraft 
entering  service. 

2. PRODUCTION. TIME  - This  real  variable  defines  the  second  aircraft 
production  rate.  It  is  the  time  in  flight  hours  between  aircraft 
entering  service. 

PRCHG  - This  real  variable  is  the  simulation  time  when  the  second 
aircraft  production  rate  takes  effect.  Note  that  this  time  is 
measured  from  the  time  that  the  first  aircraft  enters  service  and 
not  from  the  start  of  the  simulation. 

START. TEST  - This  real  variable  is  the  time  in  flight  hours  relative 
to  the  start  of  the  simulation  when  the  fatigue  test  of  an  element 
is  begun.  If  no  fatigue  test  is  to  be  conducted,  this  variable  is 
set  to  the  machine  upper  limit. 

TEST. ACCEL. FACT  - This  real  variable  is  the  fatigue  test  acceleration 
factor,  that  1 s’,  the  quotient  of  the  equivalent  flight  hours  divided 
by  the  fatigue  test  hours. 

C. GROWTH. RATE  - This  real  variable  is  the  corrosion  area  growth  rate 
in  square  inches  per  hour  for  the  aircraft  being  considered.  The 
growth  rate  for  each  element  in  the  aircraft  is  modified  by  its 
associated  CRR  (corrosion  resistance  rating). 
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C7  - If  a modification  is  developed  because  of  a fatigue  test 
failure,  this  real  variable  Is  the  percentage  (expressed  as  a 
decimal  fraction)  of  trie  test  life  when  the  inspection  frequency 
is  increased. 

C28  - This  real  variable  is  the  percentage  (expressed  as  a 
HeFimal  fraction)  reduction  in  the  remaining  fatigue  life  of  an 
element  when  corrosion  occurs  in  a stress  concentration. 

C29  - This  real  variable  is  the  percentage  (expressed  as  a 
cfeciinal  fraction)  reduction  in  the  remaining  fatigue  life  of 
an  element  when  corrosion  occurs  outside  a stress  concentration. 

MU.R  - This  real  variable  is  the  mean  of  the  log-normal 
distribution  of  the  ratio  of  the  actual  average  fatigue  life 
to  the  predicted  average  fatigue  life. 

SIG.R  - This  real  variable  is  the  standard  deviation  of  the 
log-normal  distribution  of  the  ratio  of  the  actual  average 
fatigue  life  to  the  predicted  average  fatigue  life. 

DLL  - This  real  variable  is  the  design  limit  load  in  g's  above 
the  1-g  level . 

1ABCD(1)  - This  real  variable  is  the  inspection  interval  in 
flight  hours  of  the  A-level  inspection.  It  remains  constant 
throughout  the  simulation. 

1 ABCD(2)  - This  real  variable  is  the  inspection  interval  in 
Flight  hours  of  the  B-level  Inspection.  It  remains  constant 
throughout  the  simulation. 

CABCD(*)  - This  one-dimensional  real  array  of  size  four  contains 
"tile  inspection  cost  at  each  level  of  Inspection.  CABCD(l) 
corresponds  to  the  A-level  cost;  CABCD(2)  corresponds  t.o  the 
B-level  cost;  CABCD(3)  corresponds  to  the  C-level  cost;  and 
CABCD(4)  corresponds  to  the  D-level  cost. 

S.OPT  - This  alpha  variable  is  "YES"  if  the  random  number  seeds 
are  to  be  input;  it  is  "NO"  if  seeds  are  not  input. 

LONG. LIST  - This  alpha  variable  is  "NO"  if  standard  output  and 
,rYES"  for  detailed  output  on  particular  aircraft.  (LONG  LIST 
DATA). 

FAIL. OPT  - This  integer  variable  is  "1"  for  output  of  probability 
of  failure  that  is  based  on  averaging  individual  element  failure 
rates  to  obtain  element  type  failure  rates.  This  integer  variable 
is  "2"  for  output  of  probability  of  failure  that  is  based  on  using 
a log-normal  crack  distribution  and  a curve  fit  of  probability  of 
failure  versus  crack  length  to  obtain  element  type  failure  rates. 
This  integer  variable  is  "3"  if  both  options  are  desired. 
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FAT. TEST. FACTOR  - This  real  variable  is  the  probability  of 
a fatigue  test  being  done  on  the  structural  element.  This 
real  variable  is  compared  with  a random  number  to  determine  if 
the  fatigue  test  is  done.  If  the  fatigue  test  is  not  done, 
the  fatigue  test  life  is  set  to  9999999. 

ACTUAL. AVG. FAT. LIFE  - This  real  variable  is  the  actual  average 
fa ti gue”T T f e i n fTTght  hours  determined  by  fatigue  test.  If 
this  value  is  not  known,  input  zero  and  SAIFE  will  determine 
it  statistically. 

LEAD.TIME  - This  real  variable  is  the  time  in  flight  hours 
between  when  a decision  is  made  to  develop  a structural 
modification  and  the  time  the  modification  is  available  for 
installation. 

T. FREQ.CHG  - This  real  variable  is  the  percentage  expressed 
as  a decimal  fraction  that  the  D-level  inspection  interval  is 
reduced  due  to  a fatigue  test  failure. 

S. FREQ.CHG  - This  real  variable  is  the  percentage  expressed 
as  a decimal  fraction  that  the  D-level  inspection  Interval  is 
reduced  due  to  service  experience. 

FREQ. DECREASE  - This  real  variable  is  the  percentage  expressed 
as  a dec Ima Infract ion  that  the  C and  D-level  inspection 
intervals  are  increased  due  to  favorable  service  experience. 

A.  REPAIR. COST  - This  real  variable  is  the  repair  cost  at  the 
A-level  inspection. 

B.  REPAIR. COST  - This  real  variable  is  the  repair  cost  at  the 
B-level  inspection. 

C.  REPAIR. COST  - This  real  variable  is  the  repair  cost  at  the 
C-Tevel  inspection. 

D.  REPAIR. COST  - This  real  variable  is  the  repair  cost  at  the 
D-level  inspection. 

1ST. TOOLING  - This  real  variable  is  the  tooling  cost  of  the 
first  structural  modification. 


AD. TOOLING  - This  real  is  the  tooling  cost  in  the  development 
of  any  additional  structural  modifications. 

1ST.MD.C0ST  - This  real  variable  is  the  installation  cost  of 
the  first  structural  modification. 
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AD. MD. COS'!  - This  real  variable  is  the  installation  cost  of 
any  additional  structural  modification. 

S.RF.PAIR.COST  - This  real  variable  is  the  repair  cost  of  a 
defect  detected  during  a special  Inspection. 

SU  - This  real  variable  is  the  element  ultimate  strength  in 
i^s  above  the  1-g  level. 

SF  - This  real  variable  is  the  element  fail  safe  strength  in 
g^s  above  the  1-g  level. 

SI  - SIC  - These  integer  variables  are  the  ten  random  number 
seeds  and  are  input  only  if  S.QPT  - "YES".  Any  integer  value 
may  be  used  as  input. 

LONG  LIST  DATA.  Occasionally  in  the  standard  output,  elements  will  appear 
with  unusually  long  fatigue  cracks  or  early  element  failures.  It  is 
desirable  to  have  a more  complete  service  history  of  aircraft  with  these 
early  element  failures  than  that  offered  by  the  standard  output.  This 
service  history  is  available  through  what  is  called  the  long  list  option. 
This  output  option  is  accessed  by  reading  in  alpha  characters  "YES"  for 
the  aircraft  input  variable  LONG, LIST.  After  this  input,  the  element 
description  and  identification  numbers  of  the  aircraft  to  be  tracked  are 
read  in.  The  input  variables  for  the  long  list  option  are  listed  and 
described  below  in  the  order  in  which  they  are  read  in  by  SAIFE . 

NOE  - This  integer  variable  is  the  number  of  elements  to  be 
processed  under  the  long  list  option- 

ELID(*,*)  - This  two-dimensional  alpha  array  of  size  four  by 
NOE  identifies  each  element  to  be  processed.  This  identifica- 
tion must  appear  in  the  first  sixteen  columns  of  the  data  card 
and  must  be  identical  to  the  description  read  into  the  variable 
ELEMENT (*)  described  in  ELEMENT  DATA. 

NQAC(*)  - This  one-dimensional  integer  array  of  size  NOE  is  the 
number  of  aircraft  to  be  tracked  for  each  corresponding  element. 

TLID(*,*)  - This  two-dimensional  integer  array  of  size  NOE  by 
NOACpTcontains  the  identification  numbers  of  the  aircraft  to 
be  tracked  for  a particular  element. 

ELEMENT  DATA.  The  input  variables  which  are  unique  to  each  element  and 
must  be  read  in  for  each  element  are  listed  and  described  below  in  the 
order  in  which  they  are  read  in  by  SAIFE. 
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ELEMENTS)  - This  one-dimensional  alpha  array  of  size  four 
identifies  the  element  being  simulated.  The  total  length 
of  this  identification  cannot  exceed  sixteen  characters. 

The  SAIFE  program  distinguishes  between  elements  that  are 
pressure  loaded  and  flight  loaded.  Elements  that  are  pressure 
loaded  must  start  with  characters  FUS-  and  have  MFR-  as 
characters  five  through  eight  or  FRM-  as  characters  nine 
through  twelve.  Any  characters  can  be  used  for  the  flight 
loaded  elements. 

PREDICTED. LIFE  - This  real  variable  is  the  average  element 
fatigue  life  in  flight  hours  predicted  by  analysis.  If  the 
actual  average  fatigue  life  is  known,  this  variable  can  be 
entered  as  zero. 

Ml .MEAN  - This  real  variable  is  the  average  first  external 
crack  growth  rate  in  inches  per  flight  hour. 

M2. MEAN  - This  real  variable  is  the  average  second  external 
crack  growth  rate  in  inches  per  flight  hour. 

LGHT.TO. FAILURE  - This  real  variable  is  the  length  in  inches 
at- which  the  crack  reaches  failure  under  a 1-g  load 

CONE  - This  real  variable  is  the  crack  length  in  inches  at 
which  the  first  external  crack  growth  rate  changes  to  the 
second  external  crack  growth  rate.  (First  external  critical 
crack  length). 

FSAF. LGT  - This  real  variable  is  the  length  in  inches  at  which 
the  crack  reaches  the  fail-safe  length. 

BIRTH. DEFECT. PROBABILITY  - This  real  variable  is  the  probability 
of  a production  defect. 

CRR  - This  integer  variable  is  the  corrosion  resistance  rating. 

SDM. OCCURRENCE. RATE  - This  real  variable  is  the  occurrence 
rate  of  service  damage  per  element  per  aircraft  per  flight 
hour. 

1 .PROB  - This  real  variable  is  the  probability  of  cracks 
originating  internally. 

C.PROB  - This  real  variable  is  the  probability  of  corrosion 
originating  internally. 

INT. LVL. INSP  - This  alpha  variable  is  the  letter  identifying 
the  lowest  internal  level  inspection. 
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EXT.LVL. INSP  - This  alpha  variable  Is  the  letter  Identifying 
the  lowest  external  level  Inspection. 


MOD. TEST  - This  alpha  variable  is  input  as  "YES"  if  a 
structural  modification  Is  to  be  fatigue  tested.  Otherwise 
it  is  input  as  "NO". 

LOCATED. IN. STRESS. CON  - This  real  variable  is  the  probability 
that  there  is  corrosion  in  a stress  concentration. 

1 . CDM. OCCURRENCE. RATE  - This  real  variable  is  the  initial 
corrosion  occurrence  rate  in  occurrences  per  element  per 
aircraft  per  flight  hour. 

2.  CDM. OCCURRENCE. RATE  - This  real  variable  is  the  second 
corrosion  occurrence  rate  in  occurrences  per  element  per 
aircraft  per  flight  hour. 

CDM. RATE. CHANGE  - This  real  variable  is  the  aircraft  service 
time  in  flight  hours  when  the  second  corrosion  occurrence 
rate  takes  effect. 

L.EXT  - This  real  variable  is  the  length  in  inches  at  which 
a crack  originating  internally  becomes  external. 

M3. MEAN  - This  real  variable  is  the  average  third  external 
craclc  growth  rate  in  Inches  per  flight  hour. 

M4.MEAN  - This  real  variable  is  the  average  fourth  external 
crack  growth  rate  in  inches  per  flight  hour. 

CTWO  - This  real  variable  Is  the  crack  length  in  inches  at 
Which  the  second  external  crack  growth  rate  changes  to  the 
third  external  crack  growth  rate.  (Second  external  critical 
crack  length). 

CTHREE  - This  real  variable  is  the  length  in  inches  at  which 
the  third  external  crack  growth  rate  changes  to  the  fourth 
external  crack  growth  rate.  (Third  external  critical  crack 
length) . 

I. NT. CONE  - This  real  variable  is  the  length  in  inches  at  which 
the  first  internal  crack  growth  rate  changes  to  the  second 
internal  crack  growth  rate.  (First  internal  critical  crack 
length). 

I NT. CTWO  - This  real  variable  is  the  length  in  inches  at  which 
tKi~second  internal  crack  growth  rate  changes  to  the  third 
internal  crack  growth  rate.  (Second  internal  critical  crack 
length). 
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IN.CTHREE  - This  real  variable  Is  the  length  In  Inches  at 
which  the  third  internal  crack  growth  rate  changes  to  the 
fourth  Internal  crack  growth  rate.  (Third  Internal 
critical  crack  length). 

1ABCD(3)  - This  real  variable  is  the  Initial  Inspection 
Interval  in  flight  hours  of  the  C-level  inspection. 

1ABCD(4)  - This  real  variable  is  the  initial  inspection 
Interval  in  flight  hours  of  the  D-level  inspection. 

POP. SIZE  - This  integer  variable  Is  the  number  of  elements 
of  the' same  type  on  the  aircraft.  It  is  not  necessary  to 
input  all  elements  of  the  same  type.  A sampling  may  be 
used  and  SAIFE  will  extrapolate  the  probability  of  failure 
calculation  to  the  total  number  of  elements  actually  in  the 
aircraft. 

AMEAN  - This  real  variable  is  the  result  of  fitting  an 
exponential  curve  to  flight  or  pressure  load  exceedance 
data.  AMEAN*exp(BL)  is  the  number  of  loads  per  hour  which 
exceed  the  load  level  (L). 

3_  - This  real  variable  is  the  result  of  fitting  an  exponential 
curve  to  flight  or  pressure  load  exceedance  data.  AMEAN*exp(BL) 
is  the  number  of  loads  per  hour  which  exceed  the  load  level  (L). 


FORMAT  SPECIFICATIONS,  Most  of  the  Input  data  are  entered  Into  SAIFE  by 
the  free-form  read  statement.  The  program  has  only  three  formatted  read 
statements. 

The  aircraft  type  identification,  the  alpha  array  MODEL,  is  entered 
under  the  format  specifications  2A4.  This  identification  must  be  contained 
in  the  first  eight  columns  of  the  first  card  of  the  Aircraft  Input  Data. 

All  subsequent  data  in  this  section  can  appear  in  any  columns  and  on  as 
many  cards  as  desired.  All  input  values  must  be  separated  from  one 
another  by  at  least  one  blank  column  and  a value  cannot  be  continued  on 
the  next  card. 

The  element  identification  alpha  array  ELID  is  entered  under  the  format 
specification  4A4.  This  identification  must  be  contained  in  the  first 
sixteen  columns  of  the  long  list  element  data  card.  Subsequent  data  can 
appear  in  any  column  and  on  as  many  cards  as  desired.  When  a second 
element  is  to  be  identified,  its  description  must  again  appear  in  the 
first  sixteen  columns  of  the  data  card. 

The  element  identification  alpha  array  ELEMENT  is  entered  under  the 
format  specification  4A4.  This  identification  must  be  contained  in  the 
first  sixteen  columns  of  the  first  card  of  each  set  of  Element  Input 
Data.  As  in  the  Aircraft  Input  Data,  all  subsequent  data  in  this 
section  can  appear  in  any  columns  and  on  as  many  cards  as  desired. 

SAMPLE  INPUT.  Sample  input  data  consisting  of  aircraft  input  data  and 
four  sets  of  element  input  data  are  illustrated  in  Figure  2.  The 
aircraft  type  identification  is  HYBRID.  The  four  cards  immediately 
following  the  aircraft  type  contain  the  aircraft  input  data.  The 
element  input  data  begins  immediately  after  the  element  identification 
on  the  same  card  and  terminates  on  the  last  card  before  the  next  element 
id,  itification.  The  card  following  the  last  set  of  element  input  data 
must  contain  EOD  in  the  first  three  columns.  Note  this  input  is  for 
standard  output  since  LONG. LIST  = "NO". 

Sample  input  data  to  obtainlong  list  output  on  the  last  element  (FUS- 
MRF-SID-1 740)  of  the  sample  case  are  shown  in  Figure  3.  Long  list  data 
are  requested  on  aircrafts  148  and  162.  In  this  case  we  are  duplicating 
the  output  of  the  sample  case  in  Figure  2 by  inputting  the  random  number 
seeds  appropriate  to  element  FUS-MFR- SID- 1740.  Note  that  S.OPT  = "YES" 
and  LONG. LIST  = "YES"  in  the  aircraft  data  in  Figure  3.  The  ten  random 
number  seeds  follow  the  aircraft  data.  Next  comes  the  long  list  data 
and  finally  the  element  data.  Element  output  can  be  duplicated  on  a 
given  element  without  rerunning  all  the  elements  ahead  of  the  desired 
element  by  using  the  random  number  seed  option. 
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Figure  3.  Sample  Input  Data  for 
Long  List  Output 
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PROGRAM  OUTPUT 


Each  element  to  bo  simulated  by  SAIFE  is  identified  by  three  groups  of 
alpha  characters  and  one  group  of  numeric  characters.  The  alpha  characters 
define  the  basic  element  type  and  general  location  on  the  aircraft,  and 
the  numeric  characters  define  the  specific  location  of  the  element  by 
identifying  the  wing  or  fuselage  station  number.  For  example,  an 
element  identified  as  WNG-STR-CEN-396  would  be  a wing  stringer  located 
midway  between  the  front  and  rear  spars  and  centered  at  wing  station 
396. 

m 

The  standard  program  output  consists  of  two  parts.  The  first  part 
consists  of  the  simulation  results  for  each  specific  element.  This  part 
is  printed  for  each  set  of  element  input  data.  The  second  part  consists 
of  a summary  of  the  first  parts  for  an  element  type.  In  the  example 
discussed  above,  WNG-STR-CEN  Identifies  the  element  type.  Whenever  the 
program  encounters  a set  of  element  input  data  in  which  any  single 
character  of  the  three  groups  of  alpha  characters  differs  from  those  in 
the  previous  set  of  element  input  data,  a summary  is  printed. 

The  program  also  has  a long  list  form  of  output.  This  output  consists 
of  a detailed  time  history  of  individual  elements  and  aircraft.  This 
form  of  output  would  generally  be  used  after  a standard  output  has  been 
obtained  on  element  data.  Generally,  only  the  elements  that  long  list 
output  is  desired  on  is  rerun  using  the  seed  option  to  input  the  appropriate 
random  number  seeds  for  the  desired  elements. 

STANDARD  OUTPUT,  ELEMENT  DATA.  Figure  4 illustrates  the  standard  output 
for  the  Tnput  shown  in  Figure  2.  The  random  number  seeds  and  aircraft 
numbers  for  aircraft  having  corrosion,  production  defects,  and  service 
damage  are  printed  on  the  first  page.  The  random  number  seeds  are  the 
initial  seeds  for  each  element.  Since  the  initial  seeds  for  an  element 
are  dependent  on  the  history  of  the  preceding  element,  it  is  necessary 
to  know  the  initial  element  seeds,  if  one  wanted  to  duplicate  a particular 
element  without  rerunning  all  the  previous  elements. 

On  the  next  page,  the  aircraft  type  identification,  number  of  aircraft 
in  the  fleet,  aircraft  service  life,  structural  element  identification, 
predicted  average  fatigue  life,  actual  average  fatigue  life,  and  fatigue 
test  life  are  printed  at  the  top  of  the  page. 

The  number  of  occurrences  in  the  fleet  and  the  times  to  initiation  of 
the  four  types  of  aircraft  defects  considered  by  SAIFE  are  displayed 
next.  Whenever  there  is  a fatigue  crack  initiation  and  there  are  no 
other  cracks  in  the  element,  a first  carck  is  said  to  have  occurred.  A 
single  element  can  experience  more  than  one  first  crack  in  its  lifetime 
by  having  a crack  initiation  after  a repair.  Similarly,  there  can  be 
more  than  one  occurrence  of  corrosion  and  service  damage  defect  to  exist 
simultaneously  in  the  same  element. 
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AIRCRAFT  TYPE:  HYBRID 


CRACK  LENGTHS  AND  CORRESPONDING  CUMULATIVE  PROBABILITY  OF  FAILURE 


AIRCRAFT  AO.  FL1.  HOURS  CRK.LGT.  PROB.  OF  FAILURE 


I > I t I I I I I I I I I I I I I I I I I I I I I I t t < I J I I I I I t I I I I I I t I I I I I • | I I 

UJUiUJUJUJUJUjUJUJUJUJUJkJUJa.'U.  UJUiUJUJLJLJUJUJU,UJUJUJUJlL>LJUJUjUJUJUjUJUJUJUJUJUJUJUjUjLUUJbJUJUiUjlUkJ 


,-,KnNn^'i)-#mocu,ir-.c*^-<rorvj£?>«o«j'0'~'(TrnL/t«aNUi<u--^iro^oaDOf\Jo»Nr)n^nMoOK®ao(r^ooc\j 

Nir(vo9ooinNo«9(vs<o»*A^0«<tivM/innKOf)M*i4«H4«(M*)nri9>tfi«K4’in^h>oo«(VH 


fV  <\j  m ro  -<  m # oj 


rvi  H*nj  ♦ «-4  «-<  n -t 


nj»or)Ho«oooooooooooooeciooooaoooceR)aaoooooeooaoooooooncD 

«ooHiru)invn/iiAOifieooinooeiAet/ttfiiniAiririfUflioiriAnoo«oooooooooooo0ooi«in 

oonio^(oiviNM(V4'*oeo««o««tniM«U)«NNr'9hON«i/iooooeooeeoooooaooeo«H 

«4<MA4M>*NN«««<fliAtniA«4niviNoi'>«e(Mi«4B»eDio«oo«ooooooeooeoeeeo«in 


— o n ->  -a  m ■*  ■*  * MH(yj(M(\jH'#ifl^N(o©o,wNMNh‘p<W(,itflo  -o' NN  rtnwnf 

mi  mi  mi  ^ mi  mi  mimirnir-imimiGl&ltVMt^  mi  mi 


5«3a  1.69  1.1E-13 


1 


l 


7?7777T777T7777T7777777T77777??7777777777  t 7 

UJlUUiUJUJUJUJU.kJUJUJkJMJU.UiU.LLlijU.Uj  ki  U.  u.  u.  UJ  lu  k U.  li.  Uj  U.  UjUUJU  La.U.UjUJujUjUjuJ 


fy  w m ^ « o 


H H n H 


oru>*(\jAji^einoivt«« 


rt  iy  N N 


x 

UJ 


Uj  O I 

5SX  I 


V»  CL  i 
*-*  I 
-J  4C  I 


Production  defects  are  one-time  occurrences  unless  there  is  a structural 
modification  installed.  These,  too,  can  have  production  defects.  The 
times  to  defect  initiation  are  measured  from  the  time  when  the  aircraft 
enters  service  for  the  initial  defects  and  from  the  time  when  the 
aircraft  was  last  repaired  for  subsequent  defects. 

Next,  the  number  and  lengths  of  cracks  detected  in  the  fleet  at  each  level 
of  inspection  are  printed.  These  numbers  include  second  crack  detections. 
Following  the  crack  detection  output  are  the  number  and  areas  of  corrosion 
defects  detected  in  the  fleet  at  each  level  in  inspection.  A history 
of  the  inspection  interval  changes  is  printed  next.  Each  time  that  the 
aircraft  service  experience  indicates  that  either  an  interval  increase 
or  an  interval  decrease  is  needed,  the  new  interval  values  are  printed. 

The  history  of  the  inspection  intervals  includes  the  modification  in 
effect  at  a given  Interval,  the  D-level  internal  sampling  change,  and  the 
time  the  interval  changed.  A modification  number  of  0 Indicates  the  element 
was  not  modified.  Although  the  number  of  interval  changes  allowed  in  the 
simulation  is  unlimited,  the  output  array  size  limits  the  number  printed  to 
30. 

A table  of  crack  lengths  versus  cumulative  probability  of  failure  is 
printed  next.  This  table  includes  aircraft  number,  flight  hours,  crack 
length,  and  probability  of  failure.  Flight  hours  equal  to  the  service 
time  indicate  that  the  crack  was  never  detected  or  repaired.  Times  other 
than  this  indicate  when  the  crack  was  repaired  or  a modification  was 
installed. 

The  number  of  fleet-wide  special  inspections  performed  is  printed  next. 

Next,  the  number  of  structural  modifications  developed  is  printed.  This 
number  includes  modifications  because  of  fatigue  test  failures  or  aircraft 
service  experience.  The  final  actual  average  fatigue  life  is  printed  next. 
If  there  have  been  no  modifications,  this  number  will  be  the  same  as  that  at 
the  top  of  the  page.  If  there  have  been  modifications,  this  number  is  the 
actual  average  fatigue  life  of  the  most  recently  developed  modification. 

Shown  next  is  the  number  of  aircraft  modified  in  service.  If  the  only 
modification  developed  was  due  to  a fatigue  test  failure,  this  number  can 
be  zero  if  the  test  life  was  such  that  retrofits  were  not  required.  If 
there  were  more  than  one  modification  requiring  retrofits,  this  number  can 
be  greater  than  the  size  of  the  fleet.  Printed  next  is  the  element  failure 
rate.  Finally,  each  time  an  aircraft  experiences  structural  failure  or 
its  residual  strength  reaches  its  fail-safe  strength,  the  aircraft  number 
and  the  number  of  accrued  flight  hours  are  printed.  The  aircraft  number 
is  assigned  by  its  relative  time  of  entry  into  service.  Aircraft  No.  1 
is  the  first  aircraft  to  enter  service. 
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STANDARD  OUTPUT  SUMMARY  DATA.  A suimiary  is  shown  in  Figure  4 after 
the  element  output.  Number  and  time  to  Initiation  of  aircraft 
defects,  number  and  length  of  cracks  detected  at  each  level  of 
inspection,  number  and  area  of  corrosion  defects  detected  at  each 
level  of  inspection,  initial,  mininum  and  maximum  Inspection  intervals, 
number  of  special  inspection  conducted,  number  of  structural  modifica- 
tions, number  of  aircraft  modified  in  service,  structural  failures, 
and  fail  safe  lengths  are  summaries  of  the  element  data  output.  Also 
included  is  the  element  type  failure  rate.  Element  type  failure  rate 
can  be  calculated  by  two  methods.  One  is  by  simply  taking  the  average 
of  the  individual  element  failure  rates  and  multiplying  by  the  total 
element  type  POP. SIZE.  The  other  method  uses  a log  normal  distribution 
to  generate  additional  cracks  to  account  for  the  POP. SIZE.  The 
probability  of  failure  from  the  additional  cracks  is  obtained  from  a 
semi-log  curve  fit  of  crack  length  versus  probability  of  failure  for 
the  existing  input  elements  of  the  element  type.  The  constants  of  the 
curve  fit  and  the  mean  and  standard  deviation  of  the  crack  lengths 
are  shown.  The  output  of  the  type  of  probability  of  failure 
calculation  is  dependent  on  the  input  quantity  FAIL. OPT. 

After  the  last  element  summary,  the  aircraft  failure  rate  is  printed. 
The  aircraft  failure  rate  is  based  on  pressure  and  flight  loaded 
structure.  Two  types  of  failure  rates  can  be  printed,  depending  on 
the  input  quantity  FAIL. OPT.  Next  the  five  maximum  crack  lengths 
divided  by  the  fail-safe  length  are  printed  for  flight  and  pressure 
loaded  structure  and  the  two  failure  options.  The  cracks  based  on 
averaging  are  weighted  by  the  POP. SIZE  divided  by  the  total  number 
of  element  types  that  are  input.  Thus,  the  same  crack  may  appear  more 
than  once  for  the  averaging  option.  Cracks  based  on  the  curve  fit 
option  include  both  the  cracks  that  occur  in  the  structural  element 
and  additional  cracks  that  are  generated  to  account  for  POP. SIZE. 
Finally,  a failure  summary  is  printed. 

LONG  LIST  DATA.  If  the  long  list  option  is  used,  a detailed  time 
history  of  elements  and  aircraft  are  printed.  A long  list  printout 
for  the  input  in  Figure  3 is  shown  in  Figure  5.  Long  list  print 
information  is  contained  in  the  following  events  and  routines: 


1. 

Event 

ENTER. SERVICE 

2. 

Routine 

INSTALL. MODIFICATION 

3. 

Event 

IN. SERVICE. DAMAGE 

4. 

Event 

CORROSION 

5. 

Event 

1. STRENGTH. REDUCTION 

fi. 

Event 

2. STRENGTH. REDUCTION 

7. 

Event 

l.ITE 

8. 

Event 

2.ITE 

9. 

Etent 

D. LEVEL. INSPECTION 

10. 

Routine 

EXAMINE 

11. 

Event 

REACH. FAIL. SAFE. LGT 

12. 

Event 

FAILURE 

13. 

Event 

RETIRE. FROM. SERVICE 

14. 

Event 

REPAIR 

15. 

Event 

T. INSPECTION. INCREASE 

16. 

Event 

INCREASE. INSPECTION. FREQUENCY 

17. 

Event 

IMMEDIATE. FLEET. INSPECTION 

18. 

Event 

DECISION. ON. MOD 

19. 

Event 

IMPLEMENT. MODIFICATION 
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APPENDIX  A 

DETAILED  PROGRAM  DESCRIPTION 


In  the  following  detailed  description  of  the  SAIFE  program,  each  event 
and  routine  in  the  program  is  presented  separately.  Each  presentation 
consists  of  a description,  the  definition  of  the  local  variables,  if 
any,  and  major  logic  steps  including  subroutine  calls  and  event  scheduling 
as  necessary.  All  events,  subroutines,  and  variables  are  capitalized. 
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PREAMBLE 


DESCRIPTION  - The  PREAMBLE  is  the  definition  section  of  a SIMSCRIPT 
program.  /Til  global  variables  and  global  arrays  are  defined.  Temporary 
entities  are  defined  and  tally  statistics  are  identified.  Event 
notices  and  functions  are  defined  and  an  event  priority  order  is  set. 
Global  variables  which  are  input  variables  are  not  included  here  but 
can  be  found  in  the  input  section. 

Global  Real  Variables 

FLEET. STR. RED  - This  variable  is  the  sum  of. crack  lengths  found  in 
the  fleet  since  the  last  inspection  frequency  change. 

1AAFL  - This  variable  is  the  actual  average  fatigue  life  of  the 
element  design  determined  in  the  MAIN  program. 

CRRF  - Assigned  a value  in  routine  INITIALIZATION  according  to 
the  element  corrosion  resistance  rating,  this  variable  is 
multiplied  by  the  aircraft  corrosion  growth  rate  to  give  the  element 
corrosion  growth  rate. 

COST. OF. REPAIRS  - This  variable  is  the  sum  of  repair  costs  for  the 
fleet  since  the  last  modification. 

FIX  IT. COST  - This  variable  is  the  cost  of  repairing  a defect  found 
at  a particular  inspection  level.  Its  value  is  set  in  the  inspection 
events. 

CHG. FREQ. TIME  - This  variable  is  set  equal  to  TIME.V  whenever  an 
inspection  interval  change  is  scheduled  in  the  event  REPAIR. 

1CRKT  - Each  time  a first  crack  occurs,  this  variable  is  set  equal 
tcTtlhe  service  time  on  the  aircraft. 

1C0RT  - Each  time  corrosion  occurs,  this  variabte  is  set  equal  to 
the  service  time  on  the  aircraft. 

1 SDT  ■ Each  time  service  damage  occurs,  this  variable  is  set  equal 
to  the  service  time  on  the  aircraft. 

ACRKL,  BCRKL,  CCRKL,  DCRKL,  SCRKL  - Each  time  a crack  is  found  during 
an  A- level  , B-level,  C-level,  D-Tevel , or  Special  inspection, 
the  corresponding  variable  is  set  equal  to  the  crack  length. 

ACA,  BCA,  CCA,  DCA,  SCA  - Each  time  corrosion  is  found  during  an 
A-level,  B-level,  C-level,  D-level,  or  Special  inspection,  the 
corresponding  variable  is  set  equal  to  the  corrosion  area. 


- 36  - 


variable  Is  the  number  of  flight  hours  accumulated 
since  the  last  modification  for  aircraft  no  longer  in  service, 

f^fce't1meho5iUeaatjcrraf"aCk  °CCUr5•  th<S  Variab,e  is  Set  equa’  *°  the 

servTce^time^on^ the^lrcraf t . occurs ’ th,s  1*  »«  «*»1  to  the 

ttelJrerJr9*  °CCUrS'  tMS  Variab,°  ’'S  SCt  eq“al  t0  the 

tS~ R-?p5pT~rf * G?CRn,iGSCy-  " Each  time  a crack  is  found  during  an  A- 

level  .B-level^^  evel , D-level,  or  Special  inspection,  the  correspondinq 
variable  1S  set  equal  to  the  crack  length.  currcsponoing 

GgC^_GCCA,  GDCA,  GSCA  - Each  time  corrosion  is  found  during  an  A- 

!? 5 J,  - If!  • „C-'T  • Jrltvel,  or  special  inspection,  the  ?o"e^d?pg 
variable  is  set  equal  to  the  corrosion  area, 

~S.k»  PINSL  - Each  time  there  is  an  inspection  interval  change,  these 
variables  are  set  equal  to  the  C-level  and  D-level  intervals,  respectively. 

ABC^OLD  - This  variable  is  used  to  store  the  value  of  the  C-level  inspection 
interval  prior  to  a interval  reduction.  P 

ABCD.OLD. D - This  variable  is  used  to  store  the  value  of  the  D-level 
inspection  interval  prior  to  an  interval  reduction. 

SAMP. SIZE  - This  variable  is  the  number  of  elements  of  an  element  type  that 
is  i nput • 

aTrcraftTfail ure^rate. 1S  ^ t0  St°r6  intemediate  data  in  calculating  the 
ELTYP. FAIL. RATE  - This  variable  is  the  element  type  failure  rate, 
foacls " This  vanat)le  1S  the  aircraft  failure  rate  due  to  flight 

Toads ?L'RATE  ' ™S  variable  is  the  aircraft  failure  rate  due  to  pressure 

SDT_F  - This  variable  is  used  to  store  element  residual  strength. 

RLGT  -This  variable  is  used  to  store  crack  lengths  in  calculating 

element  residual  strength  reduction  rates. 

AEXIP  - This  variable  is  the  result  of  fitting  a curve  AEXP*EXP(BA*CL)  to 
crack  length  versus  probability  of  failure. 


BA  - This  variable  Is  the  result  of  fitting  a curve  AEXP*EXP(BA*CL) 
to  the  crack  length  versus  probability  of  failure. 

AVGL  - This  variable  Is  the  average  of  element  type  crack  lengths. 

STDL  - This  variable  Is  the  standard  deviation  of  element  type 
crack  lengths. 

The  global  integer  variables  are  listed  next.  Again,  input 
variables  are  not  included  in  this  list. 

Global  Integer  Variables 

_ID  - In  each  event  and  routine,  this  variable  is  the  identification 
number  of  the  aircraft  being  processed. 

IDCK  - This  variable  is  initialized  to  zero  and  incremented  by 
one  each  time  an  aircraft  enters  service. 

X - This  variable  is  used  as  a local  index  or  array  subscript  in 
different  locations  in  the  program. 

COUNT. ELEMENT  - Each  time  new  element  data  is  read  in,  this  variable 
is  incremented  by  one. 

NICHG  - This  variable  is  the  number  of  times  that  the  inspection 
intervals  have  changed. 

LHTA  - This  variable  is  the  identification  number  of  the  aircraft 
among  the  ten  high-time  aircraft  with  the  fewest  flight  hours. 

1 . NUM. OF. RETIRE  - This  variable  is  the  number  of  aircraft  that 
have'  been  retired  from  service. 

2.  NUM. OF .CRASH  - This  variable  is  the  number  of  aircraft  which 
have  been  removed  from  service  because  of  structural  failure. 

ITRNL  - This  variable  is  the  numeric  identification  of  the  lowest 
internal  level  of  inspection. 

EXT. INSP. LEVEL  - This  variable  is  the  numeric  identification  of 
the  lowest  external  level  of  inspection. 

LIL,  LEL  - These  variables  are  the  numeric  identifications  of  the 
lowest  internal  and  external  levels  of  inspection,  respectively. 

If  either  of  these  variables  is  less  than  three,  it  is  set  equal 
to  three. 

TQ.BE.MODIF'IED  - This  variable  is  the  number  of  aircraft  with 
a pending  retrofit  modification. 
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BEEN. MODI F IED  - This  variable  is  the  number  of  aircraft  that  have 
had  a current  retrofit  modification  Installed. 


FDCK  - This  variable  is  the  number  of  aircraft  in  service  when 
a modification  is  implemented  because  of  a fatigue  test  failure. 

01 ICR.  OCOR,  OSDM,  OPD  - These  variables  are  the  number  of  occurrences 
of  first  cracks,  corrosion,  service  damage,  and  production  defects, 
respectively,  for  a particular  element. 

OSCR,  OSCO  - These  variables  are  the  number  of  cracks  and  corrosion 
defect s , Te spectively,  detected  during  a special  inspection  for  a 
particular  element. 

NSIC  - This  variable  is  the  number  of  special  Inspections  conducted 
for  a particular  element. 

NSMD  - This  variable  is  the  number  of  aircraft  modified  in  service 
for  a particular  element. 

NSFL  - This  variable  is  the  number  of  aircraft  experiencing 
structural  failure  for  a particular  element. 

NMD  - This  variable  is  the  number  of  structural  modifications 
macTe  on  a particular  element. 

NRFS  - This  variable  is  the  number  of  aircraft  with  the  residual 
strength  for  a particular  element  reaching  the  fail-safe  strength. 

SNRFS  - This  variable  is  the  number  of  aircraft  with  the  residual 
strength  for  a particular  element  type  reaching  the  fail-safe 
strength. 

-•  This  variable  is  used  as  a local  index  or  array  subscript  in 
different  locations  in  the  program. 

LDX  - If  the  long  list  option  is  in  effect,  this  variable  is  the 
ascending  sequential  position  of  the  element  being  processed 
among  those  elements  read  In  under  the  long  list  option. 

G01CR,  GOCR,  GOSDM  - These  variables  are  the  number  of  occurrences 
of  first  cracks,  corrosion  damage,  and  service  damage, 
respectively,  for  a particular  element  type. 

GOSCR , GOSCO  - These  variables  are  the  number  of  cracks  and  corrosion 
defects,  res'pectively,  detected  during  a special  Inspection  for  a 
particular  element  type. 

GOPD  - This  variable  is  the  number  of  occurrences  of  production 
defects  for  a particular  element  type. 


39  - 


SNSIC  - This  variable  is  the  number  of  special  Inspections  conducted  for 
a particular  element  type. 

SNMP  - This  variable  is  the  number  of  structural  modifications  made  on  a 
particular  element  type. 

SNSMD  - This  variable  is  the  number  of  aircraft  modified  in  service  for  a 
particular  element  type. 

SNSFL  - This  variable  is  the  number  of  aircraft  experiencing  structural 
failure  for  a particular  element  type. 

CR.CTR  - This  variable  counts  the  number  of  detected  cracks. 

OLD. SAMP  - This  variable  stores  the  value  of  D-level  internal  sampling 
pri or  "to  an  interval  reduction. 

MOD. NO  - This  variable  counts  the  number  of  structural  modifications  for 
each  element. 

NCZ  - This  variable  is  used  as  a counter  in  calculating  the  probability 
of " failure. 

NCZO  - This  variable  is  used  as  a counter  in  calculating  the  probability 
of  failure. 

NTIME  - This  variable  is  the  number  of  times  an  inspection  interval  change 
occurs. 

SAMPLING  - This  variable  is  the  sampling  at  the  internal  D-level  inspection. 
For  a sampling  of  4,  aircrafts  1,  5,  9 ...  receive  internal  inspections  on 
their  first  D-level  inspection.  Aircrafts  2,  6,  10...  receive  internal 
inspections  on  their  second  D-level  inspection.  Sampling  is  initialized 
to  48000  divided  by  1ABCD<4). 


Alpha  Variables 

TES. FAILURE  - This  variable  is  set  to  YES  if  there  is  a fatigue  test 
failure. 

PREVIOUSLY. MODIFIED  - This  variable  Is  set  to  YES  if  an  element  has  been 
modified. 

IMP .MOD. SCH  - This  variable  is  set  to  YES  if  the  event  IMPLEMENT. MODIFICATION 
is  scheduled. 

LTHO  - This  variable  is  set  to  YES  if  LONG. LIST  is  input  as  YES. 

DEC. ON. MOD. SCH  - This  variable  is  set  to  YES  if  event  DECISION. ON. MOD  is 
scheduled. 
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SELT,  SEL2,  SEL3  - These  variables  contain  parts  of  the  element  name  and 
are  used  to  distinguish  between  pressure  loaded  and  flight  loaded  structures. 

SM.CRK  - This  variable  is  set  equal  to  NO  If  a large  crack  occurs  and  YES 
if  a small  crack  occurs. 

PRIOR. CR  - This  variable  is  set  equal  to  YES  when  more  than  one  crack 
occurs  in  an  element. 

DEC. INT  - This  variable  is  set  equal  to  YES  when  an  Interval  decrease  Is 
scheduled. 

D. INT. FIND  - This  variable  is  set  equal  to  YES  if  a crack  is  found  in  an 
interval  D-level  inspection. 

DONE  - This  variable  is  set  equal  to  YES  when  the  inspection  interval  is 
first  reduced. 

IF1AG  - This  variable  is  set  equal  to  YES  when  the  inspection  interval 
Ts  reduced  due  to  a fatigue  test  failure. 


Real  Arrays 

C. INTERVAL,  D. INTERVAL  - The  elements  of  these  arrays  are  the  current 
E-TeveT  and  D-level  inspection  Intervals  for  each  aircraft  in  the  fleet. 

ABCD  - This  array  is  of  size  four  and  contains  the  most  recent  intervals 
for  each  of  the  four  levels  of  inspection. 

CKREP. TIME  - This  array  is  the  simulation  time  of  the  most  recent  crack 
repair  for  each  aircraft. 

COREP. TIME  - This  array  is  the  simulation  time  of  the  most  recent  corrosion 
repair  for  each  aircraft. 

LAST.SD  - This  array  is  the  simulation  time  of  the  most  recent  occurrence 
of  service  damage  for  each  aircraft. 

OCCUR.MOD  - This  array  Is  the  simulation  time  when  the  most  recent 
modification  was  installed  for  each  aircraft. 

SC,  SD  - These  arrays  contain  each  of  the  inspection  interval  changes  for 
the  C-level  and  D-level,  respectively. 

XRN  - This  array  is  the  random  number  selected  to  calculate  the  time  until 
structural  failure  for  each  aircraft. 

MROD  - This  array  is  the  simulation  time  of  the  most  recently  detected 
defect  at  all  inspection  levels  for  each  of  the  ten-high  aircraft. 

CGRI  - This  array  is  the  corrosion  multiplying  factor  for  each  aircraft. 


MSR1 , MFR1 , MSR2,  MFR2,  MSR1.INT,  MFR1 . INT>  MSR2.INT,  MFR2.INT  - These 
variables  are  the  Individual  element  crack  growth  rates.  There  are  four 
external  rates  and  four  Internal  rates. 

INT.LGT  - Elements  of  this  array  are  set  equal  to  L.EXT  for  aircraft 
whose  cracks  originate  Internally.  They  remain  zero  for  aircraft  whose 
cracks  originate  externally. 

SRRATE  - This  array  contains  the  aircraft  strength  reduction  rates. 

T.LAST. D - This  array  contains  the  last  D-level  inspection  Interval 
sTmuTat i on  time. 

TRCHG  - This  array  contains  times  used  in  the  calculation  of  probability 
of  failure. 

CLGT  - This  array  contains  the  length  of  all  element  type  cracks. 

AFACT  - This  array  contains  a loading  spectrum  value  drawn  from  a log 
norma!  distribution  of  AMEAN  and  divided  by  AMEAN  for  each  aircraft. 

CHG.TIME  - This  array  records  the  times  changes  are  made  in  the 
inspe-tion  intervals. 

FSAVE  - This  array  's  used  to  save  the  five  largest  cracks  on  flight 
loaded  structure  based  on  failure  option  2. 

PSAVE  - This  array  is  used  to  save  the  five  largest  cracks  on  pressure 
loaded  structure  based  on  failure  option  2. 

AVFSA  - This  array  is  used  to  save  the  five  largest  cracks  on  flight 
loaded  structure  based  on  failure  option  1. 

AVPSA  - This  array  is  used  to  save  the  five  largest  cracks  on  pressure 
loaded  structure  based  on  failure  option!. 

The  following  are  the  integer  arrays.  Again,  unless  otherwise  noted, 
all  arrays  are  1-dimensional  and  input  arrays  are  not  included. 

Integer  Arrays 

A! SR,  A2SR  - These  arrays  contain  the  event  notice  identification  numbers 
for  each  aircraft  for  the  events  1 .STRENGTH. REDUCTION,  and  2. STRENGTH. 
REDUCTION,  respectively. 

AJF  - This  array  is  the  event  notice  identification  number  for  each  aircraft 
for  event  FAILURE. 

AIRPLANE  - This  array  is  the  temporary  entity  identification  number  for 
each  aircraft. 


ML,  ABL,  ACL,  APL  - These  arrays  are  the  event  notice  identification 
numbers  f~r  each  aircraft  for  the  events  A. LEVEL. INSPECTION,  B. LEVEL. 
INSPECTION,  C. LEVEL. INSPECTION,  and  D. LEVEL. INSPECTION,  respectl vely. 

AC,  ATI  I - These  arrays  are  the  event  notice  identification  numbers  for 
each  aircraft  for  events  C0R0SI0N  and  T. INSPECTION. INCREASE,  respectively. 

ACID  - This  array  contains  the  identification  numbers  of  those  aircraft 
experiencing  structural  failure  for  a particular  element. 

OICR,  01  CO  - These  arrays  are  the  number  of  cracks  and  corrosion  defects, 
respectively,  detected  at  each  of  the  four  levels  of  inspection  for  a 
particular  element. 

SAC  ID  - This  array  contains  the  identification  numbers  of  those  aircraft 
experiencing  structural  failure  for  a particular  element  type. 

GOICR,  GOICO  - These  arrays  are  the  number  of  cracks  and  corrosion 
defects,  respectively,  detected  at  each  of  the  four  levels  of  Inspection 
for  a particular  element  type. 

HI .TIME. ACRFT  - This  array  contains  the  identification  numbers  of  the 
ten  higF-tlme  aircraft. 

API D - This  array  contains  the  identification  numbers  of  those  aircraft 
wi th  a particular  element  whose  residual  strength  has  reached  the  fail- 
safe strength. 

SAPID  - This  array  contains  the  identification  numbers  of  those  aircraft 
with  a particular  element  type  whose  residual  strength  has  reached  the 
fail-safe  strength. 

STIM  - This  array  contains  the  flight  hours  on  each  aircraft  when  the 
residual  strength  for  a particular  element  reaches  the  fail-safe  strength. 

SSTIM  - This  array  contains  the  flight  hours  on  each  aircraft  when  the 
residual  strength  for  a particular  element  type  reaches  the  fail-safe 
strength. 

FLTHR  - This  array  contains  the  flight  hours  on  each  aircraft  when 
structural  failure  occurs  for  a particular  element. 

SFLTHR  - This  array  contains  the  flight  hours  on  each  aircraft  when 
structural  failure  occurs  for  a particular  element  type. 

ARFSL  - This  array  is  the  event  notice  identification  number' for  each 
aircraft  for  event  REACH. FAIL. SAFE. LGT. 

A1E,  A2E  - These  arrays  are  the  event  notice  identification  numbers  for 
each  aircraft  for  events  l.ITE,  2.ITE,  respectively. 
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D.IN  - This  array  is  used  in  determining  which  aircraft  receive  an 
internal  D-level  Inspection. 

SAMP  - This  array  contains  the  sampling  used  at  different  inspection 
intervals. 

MOD.CTR  - This  array  counts  the  number  of  modifications  for  each 
aircraft. 

MOD. SAVE  - This  array  contains  the  modification  at  each  different 
Inspection  interval. 

PFI P - This  array  contains  the  aircraft  number  for  aircraft  that  have 
cracks. 

PFTIM  - This  array  contains  the  time  on  the  aircraft  for  cracks  used 
to  calculate  probability  of  failure. 


The  global  alpha  arrays  are  listed  next.  As  before.  Input  arrays  are 
not  Included  in  this  list. 

Global  Alpha  Arrays 

1.CR. EXISTS,  2. CR. EXISTS  - The  elements  of  these  arrays  are  set  equal  to 
"YES'*  for  each  aircraft  whenever  there  Is  a first  and  second  Initiation, 
respectively. 

CO. EXISTS  - This  array  is  set  equal  to  "YES"  for  each  aircraft  when  It 
has  corrosion  Initiation. 

SD.SCH  - This  array  is  set  equal  to  "YES"  for  each  aircraft  that  has 
event  IN. SERVICE. DAMAGE  scheduled. 

SSTAN  - This  array  is  the  station  number  which  identifies  each  aircraft 
experiencing  the  failure  of  a particular  element  type. 

SELNB  - This  array  is  the  station  number  which  identifies  each  aircraft 
with  a particular  element  type  whose  residual  strength  has  reached  the 
fail-safe  strength. 

AIL,  F$H  - This  array  is  set  equal  to  "YES"  fo>"  each  aircraft  when 
events  REACH. FAIL. SAFE. LGT  and  FAILURE,  respectively,  are  scheduled. 

IF! » IE2  - This  array  is  set  equal  to  "YES"  for  each  aircraft  that  has 
events  T. ITE  and  2.ITE,  respectively,  scheduled. 

TMOD. PENDING  - This  array  is  set  equal  to  "YES"  for  each  aircraft  that 
Fa”s"~a  modification  pending  because  of  a fatigue  test  failure. 
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SMOD . PENDING  - This  array  Is  set  equal  to  "YES"  for  each  aircraft  that 
Has  a modification  pending  because  of  service  experience. 

INSP.SCH  - This  array  Is  set  equal  to  "YES"  for  each  aircraft  that  has 
Inspections  below  the  overhaul  level  scheduled. 

l.INT,  2. I NT  - These  arrays  are  set  equal  to  "YES"  for  each  aircraft 
thatnas  aTlrst  crack  or  second  crack,  respectively,  Initiated  Internally. 

C.  INT  - This  array  Is  set  equal  to  "YES"  for  each  aircraft  that  has 
corrosion  initiated  internally. 

D. EXT  - The  appropriate  element  of  this  array  is  set  equal  to  YES  when 
an  aircraft  receives  a D-level  external  inspection. 

El,  E2,  E3,  E4  - These  arrays  are  used  to  store  the  names  of  structural 
elements  that  have  aircraft  failures. 

The  temporary  entity  definitions  and  tally  statements  are  self-explanatory. 
The  events,  functions,  and  routines  are  described  In  detail  In  the 
following  sections. 


MAIN 


DESCRIPTION  - In  the  MAIN  program,  space  is  reserved  for  all  global  arrays 
The  following  operations  are  performed  in  the  order  given:  all  Input  data 
are  read  in;  the  actual  average  fatigue  life  of  the  element  is  calculated; 
the  fatigue  test  life  is  calculated;  the  necessity  of  a structural  modifi- 
cation because  of  a fatigue  test  failure  is  determined;  the  element  type 
failure  rate  is  calculated;  the  first  event  ENTER. SERVICE  is  scheduled; 
and  the  simulation  is  initiated. 

Local  Variables 


FIRST. LIFE  - This  real  variable  is  the  fatigue  test  life  in  flight  hours. 

NFTS  - This  real  variable  is  the  time  in  flight  hours  from  which  the 
second  production  rate  goes  into  effect  to  when  the  last  aircraft  enters 
service, 

PE  - This  real  variable  is  used  to  store  intermediate  values  to  calculate 
probability  of  failure. 

PROB  - This  real  variable  is  used  to  store  Intermediate  values  to 
ca'Tc’ul ate  probability  of  failure. 

SAIJL  - This  real  variable  is  the  earliest  simulation  time  at  which  a 
structural  modification  because  of  a fatigue  test  failure  is  ready  for 
installation. 

MAJOR  LOGIC  STEPS 


1.  Reserve  aircraft  arrays. 

2.  Read  input  data. 

3.  Call  SUM. INITIALIZE. 

4.  Call  INITIALIZATION. 

5.  If  new  elelemnt  type  data: 

a.  Call  ESTIMATE. FAILURE. RATE. 

b.  Call  SUMMARY. 

6.  If  actual  average  fatigue  life  is  0:  (Call  REAL. LIFE). 

7.  Calculate  fatigue  test  life  (Call  FATIGUE. LIFE. SCATTER) . 

8.  If  fatigue  test  failure  (Schedule  T. IMPLEMENT. MOD) . 

9.  Start  simulation. 

10.  Call  DISPLAY. OUTPUT. 

11.  Read  new  element  data. 

12.  If  EOJ  STOP. 
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ROUTINE  INITIALIZATION 


DESCRIPTION  - This  routine  is  called  immediately  after  reading  each  new 
sel'  of  element  Input  by  the  MAIN  program.  This  routine  changes  the  inspection 
level  codes  to  numeric  values,  sets  the  corrosion  growth  multiplying 
factor  based  on  the  corrosion  resistance  rating,  and  resets  the  tally 
counters.  It  also  Initializes  all  the  element  global  variables  which  are 
not  part  of  the  input. 


ROUTINE  SUM. INITIALIZE 


DESCRIPTION  - This  routine  Is  called  from  the  MAIN  program  each  time  a new 
element  type  is  read  In.  The  element  type  is  identified  by  the  first 
twelve  characters  of  the  element  Identification.  This  routine  Initializes 
the  global  variables  and  rests  the  tally  counters. 


ROUTINE  COMP. RISK 


DESCRIPTION  - This  routine  determines  the  probability  of  failure  due  to  a 
first  or  second  crack.  Failure  due  to  a first  crack  is  based  on  a five 
part  strength  reduction  curve  and  failure  due  to  a second  crack  Is  based 
on  a nine  part  strength  reduction  curve.  If  a second  crack  exists  the 
lengtli  relationships  of  the  first  and  second  crack  must  be  determined. 

These  are  dependent  on  the  three  critical  crack  lengths  and  the  four  crack 
growth  rates.  The  logic  for  testing  all  the  combinations  of  first  crack 
critical  lengths  to  second  crack  critical  lengths  Is  lengthy.  This 
routine  is  called  from  routine  INSTALL. MODIFICATION  and  events  CORROSION, 

2. STRENGTH. REDUCTION,  REACH. FAIL. SAFE. LGT,  FAILURE,  RETIRE. FROM. SERVICE 
and  REPAIR. 

Local  Variables 

A - This  real  variable  is  obtained  for  each  aircraft  from  a log-normal 
distribution  defined  by  AMEAN. 

AZRD  - This  real  variable  is  used  to  store  intermediate  values  in  calculat- 
ing the  probability  of  failure. 

CCL1,  CCL2,  CCL3  - These  real  variables  are  the  first,  second,  and  third 
cri ti cal  crack'  Tengths . 

CL  - This  real  variable  is  the  sum  of  the  first  and  second  crack  lengths. 

CL1 , CL2  - These  real  variables  are  the  first  and  second  crack  lengths. 

DLC1 , DL1 , DL2  - These  real  variables  are  used  to  store  intermediate 
values  in  calculating  first  and  second  crack  length. 

K14,  K15,  K16,  K17,  K.18,  K19,  K20,  K21 , K22,K23.  K24,  K25,  K26,  K27  - 
Thes’e  real  variables  are  used  to  store  intermediate  results  in  calculating 
the  probability  of  failure. 

RS  - This  real  variable  is  the  residual  strength  of  an  element. 

SR,  SRI  , SR2,  SR3,  SR4,  SR5,  5R6,  SR'/  - These  real  variables  are  strength 
reduction  rates.  ~ 

SI,  S2 , S3,  S4,  S5,  S6  - These  real  variables  are  strength  reductions  at 
tfietfiree  critical  crack  lengths  for  the  two  cracks. 

TAC  - This  real  variable  is  the  simulation  time  to  corrosion  initiation. 

TA1  - This  real  variable  is  the  simulation  time  to  the  first  crack  Initia- 
tion. 


TA2  - This  real  variable  is  the  simulation  time  to  the  second  crack  Initia- 
tion. 


T,  Tl,  Til,  T12,  T2,  T21 , T22,  T23,  T3.  T4,  T5,  T6  - These  real  variables 
are  intermediate  time  variables  used  in  calculating  the  probability  of 
failure. 


ROUTINE  ESTIMATE. FAILURE. RATE 


DESCRIPTION  - This  routine  determines  the  failure  rate  of  an  element  type. 
TTi is  routine  fits  a semi -log  curve  to  crack  length  versus  probability  of 
failure  for  all  cracks  that  occur  in  the  individual  elements  of  the 
element  type.  If  a sample  of  elements  is  used,  the  routine  determines 
the  average  number  of  cracks  per  element  and  extrapolates  the  total 
number  of  cracks  to  the  actual  element  population.  The  additional  cracks 
are  generated  with  a log-normal  distribution  using  the  mean  and  standard 
deviation  of  the  cracks  in  the  element  sample.  The  probability  of 
failure  of  the  additional  cracks  is  determined  from  the  curve  fit  of 
crack  length  versus  probability  of  failure.  The  element  type  failure 
rate  is  then  determined  and  returned  to  the  calling  program.  In  addition, 
the  five  largest  cracks  in  terms  of  fail-safe  length  are  determined  and 
saved  each  time  the  routine  is  called.  Cracks  are  saved  In  four 
categories  as  described  in  the  discussion  on  program  output.  This  routine 
is  called  from  the  MAIN  program. 

Local  Variables 


C - This  real  variable  is  used  to  count  the  number  of  cracks  with  a 
probability  of  failure  equal  to  zero. 

CL  - This  real  variable  is  the  crack  length. 

LVS,  XLYS  - These  real  variables  are  used  as  intermediate  storage  in  the 
curve  fit  procedure. 

NUM.CRKS  - This  real  variable  is  the  additional  cracks  generated  in  the 
routine  to  account  for  the  POP. SIZE  of  the  element  type. 

PE,PF,PROB,DET  - These  real  variables  are  used  as  intermediate  storage 
in  calculating  the  element  type  failure  rate. 

SSQL  - This  real  variable  is  the  sum  of  the  squares  of  the  crack  lengths. 
SUML  - This  real  variable  is  the  sum  of  the  crack  lengths. 

T - This  real  variable  is  the  fleet  service  life. 

I , IK,  J , K,  KI , KK,  NN,  TEMP  - These  variables  are  used  in  determining 
the  five  largest  cracks  to  save. 
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ROUTINE  REAL. LIFE 


DESCRIPTION  - This  routine  accepts  (1)  the  predicted  average  fatigue  life 
of  an  element  (2)  the  mean  and  standard  deviation  of  the  ratio  of  the 
actual  average  fatigue  life  to  the  predicted  average  fatigue  life.  A 
random  selection  Is  made  from  a log-normal  distribution  and  multiplied  by 
the  predicted  average  fatigue  life.  The  ratio  of  the  actual  element 
fatigue  life  to  the  predicted  fatigue  life  is  limited  from  .1  to  10.0. 

The  resulting  actual  average  fatigue  life  is  returned  to  the  calling 
program.  REAL. LIFE  can  be  called  from  the  MAIN  program  and  events 
IMPLEMENT. MODIFICATION  and  T. IMPLEMENT. MOD. 

Local  Variables 


MEAN  - This  real  variable,  whose  value  is  passed  from  the  calling  routine, 
i s the  mean  of  the  ratio  distribution. 

RATIO  - This  real  variable,  determined  to  be  log-normally  distributed,  is 
the  ratio  of  the  actual  fatigue  life  of  an  element  design  to  its  predicted 
fatigue  life. 

STD. DEV  - This  real  variable,  whose  value  is  passed  to  the  calling  routine, 
is'  the- standard  deviation  of  the  ratio  distribution. 

PPL  - This  real  variable  is  the  design  predicted  average  fatigue  life 
passed  from  the  calling  routine. 

RFL  - This  real  variable  is  the  element  actual  average  fatigue  life  which 
Ts  returned  to  the  calling  routine. 
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EVENT  ENTER. SERVICE 


DESCRIPTION  - This  event  represents  the  entry  into  service  of  a new 
aircraft.  The  temporary  entity  AIRCRAFT  is  created  and  Identified  by  the 
variable  AIRPLANE  (ID).  The  entity  attributes  TAIL. ID  and  ENTRY. TIME  are 
defined  and  the  AIRPLANE  is  filed  in  the  set  ACTIVE. FLEET.  AFACT  is 
determined  for  each  aircraft  from  a log-normal  distribution  using  input 
AMEAN.  The  routine  FATIGUE. LIFE. SCATTER  is  called  to  determine  the  times 
to  first  and  second  crack  initiation.  Element  external  and  Internal  crack 
growth  rates  are  calculated.  The  times  to  corrosion  initiation  and 
service  damage  are  calculated.  If  either  of  these  times  is  less  than  the 
service  life  of  the  aircraft,  the  corresponding  defect  is  scheduled.  The 
probability  of  a production  defect  is  determined.  If  there  Is  a pro- 
duction defect,  the  time  to  first  crack  initiation  is  replaced  by  a time 
drawn  from  a distribution  of  times  to  crack  initiation  of  aircraft  with 
production  defects.  Crack  initations,  D-level  inspection,  and  retirement 
from  service  are  also  scheduled.  If  the  present  aircraft  is  not  the  last 
aircraft  in  the  fleet,  another  ENTER. SERVICE  is  scheduled.  This  event  is 
scheduled  in  the  MAIN  program  and  within  itself. 

Local  Variables 

ASP  - This  real  variable  is  the  standard  deviation  equal  to  .15  *AMEAN 
used  to  determine  AFACT. 

DEFECT. LIFE  - This  real  variable  is  the  time  to  first  crack  initiation 
when  the  aircraft  has  a production  defect. 

FIRST. LIFE  - This  real  variable  is  the  time  to  first  crack  Initiation  when 
the  aircraft  has  no  production  defect. 

HOURS. TO. CORROSION  - This  real  variable  is  the  time  to  corrosion  initia- 
tion. 

INT. FASTI  - This  real  variable  is  the  standard  deviation  of  the  second 
Internal  crack  growth  rate  distribution. 

INT.FAST2  - This  real  variable  is  the  standard  deviation  of  the  fourth 
internal  crack  growth  rate  distribution. 

INT. SLOW!  - This  real  variable  is  the  standard  deviation  of  the  first 
internal  crack  growth  rate  distribution. 

INT. SLOW?  - This  real  variable  is  the  standard  deviation  of  the  third 
internal  crack  growth  rate  distribution. 

OURS. TO. SERVICE. DAMAGE  - This  real  variable  is  the  time  to  service  damage 
occurrence. 
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RN  - This  real  variable  is  a uniformly  distributed  random  number  between 
zero  and  one. 


STD, FASTI  - This  real  variable  is  the  standard  deviation  of  the  second 
external  crack  growth  rate  distribution. 


STD.FAST2  - This  real  variable  is  the  standard  deviation  of  the  fourth 
external  crack  growth  rate  distribution. 


STD, SLOW!  - This  real  variable  is  the  standard  deviation  of  the  first 
external  crack  growth  rate  distribution. 


STD.SL0W2  - This  real  variable  is  the  standard  deviation  of  the  third 
external- crack  growth  rate  distribution. 

MAJOR  LOGIC  STEPS 


1 . Create  AIRCRAFT 

2.  Calculate  times  to  first  and  second  crack  initiation 
(Call  FATIGUE. LIFE. SCATTER) 

3.  Calculate  external  and  internal  crack  growth  rates  (Call 
RATE  eight  times). 

4.  Calculate  corrosion  time  (Call  PREEDICT. CORROSION) . 

5.  Calculate  service  damage  time  (Call  PREEDICT. SERVICE. DAMAGE). 

6.  If  corrosion  occurs  in  aircraft  service  life  (Schedule 
CORROSION). 

7.  If  service  change  occurs  in  aircraft  service  life  (Schedule 
IN. SERVICE. DAMAGE). 

8.  Determine  if  production  defect  occurs. 

9.  If  first  crack  occurs  in  aircraft  service  life  (Schedule 
1. STRENGTH. REDUCTION). 

10.  If  second  crack  occurs  in  aircraft  service  life  (Schedule 
2. STRENGTH. REDUCTION). 

11.  Schedule  D. LEVEL. INSPECTION. 

12.  Schedule  RETIRE. FROM. SERVICE. 

13.  Schedule  ENTER. SERVICE  for  next  aircraft  if  not  last  aircraft. 


ROUTINE  FATIGUE. LIFE, SCATTER 


DESCRIPTION  - This  routine  receives  the  actual  average  fatigue  life  of  the 
element  design  and  calculates  the  fatigue  test  life  and  the  time  to  crack 
initiation  of  the  two  fatigue  cracks.  These  times  are  random  selections 
from  a two-parameter  Weibull  distribution.  If  the  routine  is  not  being 
used  to  calculate  fatigue  test  life,  then  crack  Initiation  times  are 
adjusted  by  dividing  them  by  AFACT.  This  routine  can  be  called  from  the 
MAIN  program,  INSTALL. MODIFICATION,  ENTER. SERVICE,  and  REPAIR. 

Local  Variables 


ALPHA  - This  real  variable  is  the  shape  parameter  of  the  fatigue  life 
distribution. 

BETA  - This  real  variable  is  the  scale  parameter  of  the  fatigue  life 
distribution. 

FIRST. LIFE  - This  real  variable  is  the  fatigue  test  life  or  time  to  first 
crack  initiation. 

LIFE  - This  real  array  of  length  two  is  used  to  temporarily  store  the 
times  to  crack  initiation. 

JN  - This  integer  variable,  passed  from  the  calling  routine,  identifies  the 
random  number  stream  to  be  used. 

RFL  - This  real  variable  is  the  element  actual  average  fatigue  life 
passed  from  the  calling  routine. 

RN  - This  real  variable  is  a uniformly  distributed  random  number. 

SECOND  LIFE  - This  real  variable  is  the  time  to  second  crack  intiation. 
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ROUTINE . INSTALL . MODIFICATION 


DESCRIPTION  - This  routine  represents  the  Installation  of  a structural 
modification  caused  by  a fatigue  test  failure  or  by  aircraft  service 
experience.  The  modification  Is  Installed  during  a repair  or  a D-level 
Inspection.  All  previously  scheduled  defect  Initiations  are  cancelled, 
and  new  times  to  defect  Initiations  are  calculated  for  each  aircraft  when 
it  is  modified.  This  routine  can  be  called  from  the  event  REPAIR  and 
D. LEVEL. INSPECTION. 

Local  Variables  (Same  as  ENTER. SERVICE  except  no  variable  ASD). 


MAJOR  LOGIC  STEPS 


1.  Call  COMP. RISK  If  crack  exists. 

2.  If  inspection  scheduled  (Call  CANCEL. SCHEDULED. INSPECTIONS) . 

3.  Cancel  FAILURE,  CORROSION,  1 .STRENGTH. REDUCTION,  2. STRENGTH. 
REDUCTION,  l.ITE,  2.ITE,  and  REACH. FA I L. SAFE. LGT  If 
scheduled. 

4.  Do  steps  2-10  in  ENTER. SERVICE. 
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EVENT  IN. SERVICE. DAMAGE 


DESCRIPTION  - This  event  represents  the  occurrence  of  a service  damage 
defect"  TFls  occurrence  results  in  the  immediate  Initiation  of  the  next 
scheduled  crack.  A new  time  to  service  damage  Is  determined.  If  the  new 
time  Is  less  than  the  remaining  aircraft  service  time,  this  event  Is 
scheduled  once  again.  This  event  can  be  scheduled  from  within  itself  or 
in  event  ENTER. SERVICE. 

Local  Variables 


IDSDM  - This  Integer  variable  Is  the  Identification  number  of  the  aircraft 
from  which  the  event  was  scheduled. 

OURS. TO. SERVICE. DAMAGE  - This  real  variable  is  the  time  to  the  next 
occurrence  of  service  damage. 

RST  - This  real  variable  is  the  remaining  service  time  to  the  retirement 
of  the  aircraft  being  considered. 
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EVENT  T. IMPLEMENT. MOD 


DESCRIPTION  - This  event  represents  the  development  of  a structural 
modification  because  of  a fatigue  test  failure.  A new  predicted  life  for 
the  modification  is  determined  by  using  the  old  predicted  life  or  twice 
the  actual  average  fatigue  life,  whichever  is  greater.  The  actual  average 
fatigue  life  of  the  modification  Is  determined  by  calling  REAL. LIFE.  If 
the  modification  Is  fatigue  tested  and  the  actual  average  fatigue  life  Is 
less  than  the  predicted  life,  the  actual  average  fatigue  life  Is  set  equal 
to  the  predicted  life.  If  the  fatigue  test  failure  occurs  before  one 
service  life,  a retrofit  modification  is  indicated  and  an  Inspection 
Increase  is  scheduled  for  active  aircraft.  This  event  is  scheduled  In  the 
MAIN  program. 

Local  Variables 

NEW. LIFE  - This  real  variable  is  the  predicted  life  of  the  modification. 

NMU  - This  real  variable  Is  MU.R  + .15  (1--MU.R)  ancj  based  on  the 
assumption  that  a modification  usually  improves  the  actual  average  fatigue 
life  of  a particular  design. 

NSIG  - This  real  variable  Is  .85(SIG.R)  and  is  based  on  the  same  assump- 
tion as  for  NMU. 


ROUTINE  PREEDICT. SERVICE. DAMAGE 


DESCRIPTION  - This  routine  generates  the  time  to  service  damage  occurrence 
for  a given  aircraft  from  a constant  service  damage  occurrence  rate.  If 
the  service  damage  occurrence  rate  Is  zero  in  the  Input,  the  routine 
sets  the  time  to  service  damage  to  twice  the  aircraft  service  life. 

This  routine  can  be  called  from  events  ENTER. SERVICE  and  IN. SERVICE. DAMAGE. 

Local  Variables 


OURS. TO. SERVICE. DAMAGE  - This  real  variable  is  the  time  to  service 
damage  in  flight  hours. 

RN  - This  real  variable  is  a uniformly  distributed  random  number  between 
zero  and  one. 
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ROUTINE  PREDICT. CORROSION 


DESCRIPTION  - This  routine  generates  time  to  corrosion  Initiation  for  a 
given'  aircraft  from  a time-dependent  occurrence  rate  approximated  by  two 
constant  rates.  The  first  constant  occurrence  rate,  the  second  constant 
occurrence  rate,  and  the  service  time  on  the  aircraft  when  the  second  rate 
goes  into  effect  are  all  input  variables.  This  routine  can  be  called  from 
the  routine  INSTALL. MODIFICATION  and  events  ENTER. SERVICE  and  REPAIR. 

Local  Variables 

CRCT  - This  real  variable  is  the  remaining  time  in  flight  hours  until  the 
second  corrosion  occurrence  rate  goes  into  effect.  This  variable  can  be 
negative  indicating  that  the  second  rate  is  already  in  effect. 

HOURS. TO. CORROSION  - This  real  variable  is  the  flight  time  until  corrosion 
initiation. 

LP  - This  real  variable  is  used  to  hold  an  intermediate  value  during  the 
calculation  of  time  to  corrosion  initiation. 

RN  - This  real  variable  is  a uniformly  distributed  random  number  between 
zero  and  one. 


I 


I 


i' 


i 
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EVENT  CORROSION 


DESCRIPTION  - This  event  represents  the  Initiation  of  a corrosion  defect. 

A probalistic  determination  Is  made  of  whether  corrosion  occurs  externally 
or  internally  and  If  It  occurs  at  a stress  concentration.  A corrosion 
damage  factor  is  calculated  depending  on  whether  corrosion  occurs  In  a 
stress  concentration.  A multiplying  factor  for  crack  growth  rates  Is  set 
equal  to  the  numeric  value  for  the  corrosion  resistance  rating.  If  either 
of  the  events  FAILURE  or  REACH. FAIL. SAFE. LGT  is  scheduled,  its  remaining 
time  to  occurrence  Is  reduced  by  the  corrosion  resistance  rating.  The 
remaining  time  to  crack  Initiation  is  reduced  by  the  corrosion  damage 
factor.  This  event  can  be  scheduled  in  the  routine  INSTALL. MODIFICATION 
and  in  events  ENTER. SERVICE  and  REPAIR. 

Local  Variables 

COM. MULTIPLYING  FACTOR  - This  real  variable  is  the  factor  which  when 
multiplied  by  the  remaining  time  to  crack  initiation  accounts  for  the 
shortening  effect  of  corrosion  on  fatigue  life. 

IDCO  - This  integer  variable  contains  the  identification  number  of  the 
aircraft  for  which  the  event  CORROSION  was  scheduled. 

NFTM  ~ If  a FAILURE  has  been  scheduled,  this  real  variable  is  the  remain- 
ing time  until  its  occurrence. 

REDUCED. REMAINING  LIFE  - This  real  variable  is  the  REMAINING. LIFE  multiplied 
by  the  corrosion  damage  factors. 

REMAINING. LIFE  - This  real  variable  is  the  remaining  time  until  a scheduled 
crack  initiation. 

RST  - This  real  variable  is  the  remaining  service  time  of  the  aircraft 
umTer  consideration. 

TRT  - If  a REACH. FAIL. SAFE. LGT  has  been  scheduled,  this  real  variable  is 
the  remaining  time  until  its  occurrence  multiplied  by  the  corrosion 
damage  factor. 

MAJOR  LOGIC  STEPS 

1.  Call  COMP. RISK  if  crack  exists. 

2.  Determine  if  corrosion  is  internal  or  external. 

3.  If  inspections  not  scheduled  (Call  INSPECTION. SCHEDULER) . 

4.  Determine  if  corrosion  in  stress  concentration  factor. 

5.  Calculate  corrosion  damage  factor. 

6.  Set  crack  growth  rate  multiplying  factor  to  corrosion 
resistance  rating. 

7.  Cancel  l.ITE,  2.ITE,  FAILURE,  REACH. FAIL. SAFE. LGT, 

1 STRENGTH. REDUCTION,  and  2 .STRENGTH. REDUCTION  if  scheduled. 

8.  Schedule  l.ITE,  2.ITE,  FAILURE,  REACH. FAIL-SAFE. LGT, 

1. STRENGTH. REDUCTION,  and  2. STRENGTH. REDUCTION  based  on 
time  reduction  by  corrosion  resistance  rating  or  corrosion 
damage  factor. 
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ROUTINE  RATE 


DESCRIPTION  - This  routine  statistically  generates  element  crack  growth 
rates  which  reflect  variation  in  material  properties  and  load  environment. 
The  growth  rates  are  randomly  drawn  from  a normal  distribution  which  Is 
defined  by  a mean  growth  rate  and  a standard  deviation  passed  from  the 
calling  routine.  If  a random  draw  yields  a negative  growth  rate,  the  rate 
Is  set  equal  to  the  mean  growth  rate  minus  four  standard  deviations. 

Thus,  the  user  must  be  sure  that  the  standard  deviation  is  always  less 
than  one-fourth  of  the  means.  This  routine  Is  defined  as  a function  in 
the  PREAMBLE  and  is  used  in  event  ENTER. SERVICE  and  INSTALL. MODIFICATION. 

Local  Variables 


G1  ,G2  - These  real  variables  are  used  to  hold  intermediate  values  In  the 
calculation  of  the  crack  growth  rate. 

M - This  real  variable  is  the  mean  crack  growth  rate  passed  from  the 
calling  routine. 

RN  - This  real  variable  Is  a unlformally  distributed  random  number  between 
zero  and  one. 

RNF  - This  real  variable  is  1 - RN. 

S^  - This  real  variable  is  the  crack  growth  rate  standard  deviation  passed 
from  the  calling  routine. 

W - This  real  variable  Is  used  to  hold  Intermediate  values  in  the  calcula- 
tion of  the  crack  growth  rate. 

I - This  real  variable  Is  the  element  crack  growth  rate  returned  to  the 
calling  routine. 
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EVENT  1 .STRENGTH . REDUCTION 


DESCRIPTION  - This  event  represents  the  initiation  of  the  first  crack. 

A uniform  random  number  Is  compared  with  the  probability  of  Internal 
cracking  to  determine  whether  this  crack  Initiates  Internally.  If  It 
does  Initiate  Internally,  the  time  until  It  becomes  external  is  calcu- 
lated and  the  event  l.ITE  Is  scheduled.  The  time  to  structural  failure 
is  calculated  using  a five-part  residual  strength  curve.  If  this  time 
is  less  than  the  remaining  service  life  of  the  aircraft,  the  event 
FAILURE  is  scheduled.  The  time  until  the  residual  strength  of  the 
element  reaches  the  fail  safe  length  is  calculated.  If  this  time  is 
less  than  the  remaining  service  life  of  the  aircraft,  the  event  REACH. 
FAIL. SAFE. LGT  is  scheduled.  This  event  can  be  scheduled  in  events 
ENTER. SERVICE,  INSTALL. MODIFICATION,  and  REPAIR. 

Local  Variables 


A - This  real  variable  is  obtained  for  each  aircraft  from  a log-normal 
distribution  defined  by  AMEAN. 

ARG,  Kl,  K2 , K4,  K8,  K9,  K10,  Kll , K12,  K.13,  K14,  K15,  K16,  K17,  K18, 

re , L'OKE  - THese" real  variables  are  used 
as  intermediate  1/aT ues  in  the  calculation  of  time  until  structural 
failure. 

GR1  - This  real  variable  is  the  element  first  crack  growth  rate. 

GR2  - This  real  variable  is  the  element  second  crack  growth  rate. 

GR3  - This  real  variable  is  the  element  third  crack  growth  rate. 

GR4  - This  real  variable  is  the  element  fourth  crack  growth  rate. 

ID! SR  - This  integer  variable  is  the  aircraft  identification  number. 

RJ , R2,  R3 , R4,  R5  - These  real  variables  represent  the  strength  degrada- 
tTorTon  the  five-part  strength  degradation  curve. 

SI,  S2 , S3  - These  real  variables  represent  the  residual  strengths  at 
CORETTIW,  and  CTHREE. 

TAR  - This  real  variable  is  the  simulation  time  at  which  the  aircraft 
"EeTing  processed  retires  from  service. 

T - This  real  variable  is  the  time  in  flight  hours  until  a crack  initiated 
Internally  becomes  external. 

T1 , T2 , T3,  T4  - These  real  variables  are  the  time  required  for  the 
crack"  to  grow  to  CONE,  CTWO,  CTHREE,  arid  FSAF.LGT. 

TTF  - This  real  variable  is  the  time  in  flight  hours  until  element 
fall ure. 
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MAJOR  LOGIC  STEPS 


1.  Determine  if  crack  initiates  internally. 

2.  Calculate  time  for  Internal  crack  to  become  external. 

3.  If  crack  internal  and  becomes  external  within  aircraft  service 
1 ife  (Schedule  1 . ITE) . 

4.  If  inspections  not  scheduled  (Call  INSPECTION. SCHEDULER). 

5.  Determine  time  to  reach  fail  safe  length. 

6.  If  crack  reaches  fail  safe  length  within  aircraft  service 
life  (Schedule  REACH. FAIL. SAFE. LGT) . 

7.  Determine  time  to  failure. 

8.  If  time  to  failure  within  aircraft  service  life  (Schedule 
FAILURE). 
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EVENT  2. STRENGTH. REDUCTION 


DESCRIPTION  - This  event  represents  the  second  crack  Initiation.  A 
uniform  random  number  Is  compared  with  the  probability  of  internal 
cracking  to  determine  whether  this  crack  Initiates  internally.  If  It 
does  Initiate  Internally,  the  time  until  it  becomes  external  Is  cal- 
culated and  the  event  2.ITE  Is  scheduled  If  it  occurs  within  the  service 
life  of  the  airplane.  The  time  to  structural  failure  Is  calculated 
using  a nine-part  residual  strength  curve,  which  is  dependent  on  the 
crack  growth  rate  history  of  both  the  first  and  second  crack.  If  the 
time  until  failure  is  less  than  the  remaining  service  life  of  the 
aircraft,  the  event  FAILURE  is  scheduled.  The  time  until  the  residual 
strength  of  the  element  reaches  fail  safe  strength  is  calculated.  If 
this  time  is  less  than  the  remaining  service  life  of  the  aircraft,  the 
event  REACH. FAIL. SAFE. LGT  is  scheduled.  This  event  can  be  scheduled  in 
events  ENTER. SERVICE,  INSTALL. MODIFICATION,  and  REPAIR. 

Local  Variables 


A - This  real  variable  is  obtained  for  each  aircraft  from  a log-normal 
distribution  defined  by  AMEAN. 

ARG,  Kl,  K2,  K4,  K8,  K9,  K10,  Kll , K12,  K13,  KT 4 , K1S,  Kl6,  K17,  K18, 

K19,  K20,  K21,  K22,  K23,  K24,  K2 5 , K26,  K2>,  K2¥TT29,  K30,  K31 , K32, 

K33,  K34 , K35,  K36,  oXToff,  K39 , K40,  K41 , K42,  K43,  K44,  LGK5,  LG  - 
These  local  variables  are  used  as  i ntermediate  values  in  thecal culation 
of  time  to  element  failure. 

CCL1 , CCL2,  CCL3  - These  real  variables  equal  CONE,  CTWO,  and  CTHREE 
respectively. 

DL  - This  real  variable  is  the  length  of  the  first  crack  at  corrosion 
initiation. 

GR1 , GR2,  GR3,  GR4,  GR5 , GR6,  GR7,  GR8,  GR9  - These  real  variables  are 
the  combined  growth  rates  oTthe  first  and  second  cracks. 

ID2SR  - This  integer  variable  is  the  aircraft  identification  number. 

R1 , R2,  R3,  R4,  R5,  R6,  R7,  R8,  R9,  These  real  variables  represent  the 
strength  degradation  on  the  nine-part  strength  degradation  curve. 

SI , S2,  S3,  S4,  S5,  S6,  S7  - These  real  variables  represent  the  residual 
strengths"  at  tne  actual  crack  length  and  fail  safe  length  of  the  first 
and  second  cracks. 

T - This  real  variable  is  the  time  in  flight  hours  until  a crack  initiated 
internally  becomes  external. 

TAC  - This  real  variable  is  the  time  of  corrosion  initiation. 
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TAR  - This  real  variable  is  the  simulation  time  at  which  the  aircraft 
being  processed  retires  from  service. 

TA1  - This  real  variable  is  the  time  of  the  first  crack  initiation. 

TCL1  - This  real  variable  is  the  time  for  the  crack  to  reach  CONE. 

TCLL  - This  real  variable  is  the  time  for  the  crack  to  grow  between  CONE 
and  CTWO. 

TCL3  - This  real  variable  is  the  time  for  the  crack  to  grow  between  CTWO 
ancTC  THREE. 

Tl , T2,  T3,  T4,  T5,  T6,  T7,  T8  - These  real  variables  are  the  times 
required  for  tbe  two  cracks  to  reach  critical  crack  lengths  and  fail 
safe  lengths. 

MAJOR  LOGIC  STEPS 

1.  Call  COMP. RISK. 

Determine  if  crack  initiates  internally. 

3.  Calculate  time  for  internal  crack  to  become  external. 

4.  If  crack  internal  and  becomes  external  within  aircraft 
• service  life  (Schedule  2.ITE). 

5.  Determine  time  for  two  cracks  to  reach  fail  safe  length. 

6.  If  cracks  reach  fail  safe  length  within  aircraft  service 
life  (Schedule  REACH. FAIL. SAFE. LGT) . 

7.  Determine  time  to  failure. 

0.  If  time  to  failure  within  aircraft  service  life 
(Schedule  FAILURE). 
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EVENT  1 . ITE 


DESCRIPTION  - This  event  represents  the  time  when  a first  crack  which 
Inltlatedinternally  becomes  external.  This  time  Is  defined  by  L.EXT,  an 
Input  quantity.  At  this  time  the  appropriate  element  of  the  alpha  array 
l.INT  Is  changed  from  "YES"  to  "NO".  This  event  Is  scheduled  In  event 


l.KNT  Is  changed  from  "YES"  to  "NO". 
1. STRENGTH. REDUCTION. 

Local  Variables 


ID1_E  - This  Integer  variable  Is  the  aircraft  Identification  number. 


DESCRIPTION  - This  event  represents  the  time  when  a second  crack  which 
initiated  Internally  becomes  external.  This  time  Is  defined  by  L.EXT, 
an  input  quantity.  At  this  time  the  appropriate  element  of  the  alpha 
array  2.INT  is  changed  from  "YES"  to  "NO".  This  event  Is  scheduled  in 
event  2. STRENGTH  REDUCTION. 

Local  Variables 


ID2E  - This  integer  variable  is  the  aircraft  identification  number. 


ROUTINE  INSPECTION, SCHEDULER 


DESCRIPTION  - This  routine  schedules  all  inspections  below  the  D-level  on 
a given  aircraft.  To  conserve  execution  time,  the  Inspections  are  scheduled 
so  that  the  aircraft  is  not  Inspected  before  the  defect  reaches  Its  minimum 
detectable  size  at  the  A,  B,  ana  C inspection  levels.  This  routine  is 
called  from  events  1 .STRENGTH. REDUCTION  and  C0R0SI0N. 

Local  Variables 


Cl  - This  real  variable  Is  the  corrosion  growth  rate  used  to  calculate  the 
time  to  the  minimum  detectable  corrosion  area. 

Ml  - This  real  variable  is  the  crack  growth  rate  used  to  calculate  the 
tTme  to  the  minimum  detectable  crack  length. 

N - This  Integer  variable  indicates  whether  a crack  initiation  or  a corrosion 
Initiation  caused  this  routine  to  be  called. 

S. INSP.AT  - This  real  variable  is  the  simulation  time  at  which  the  first 
Inspection  at  each  level  is  scheduled. 

TML  - This  real  variable  is  the  time  to  the  minimum  detectable  defect  size 
calculated  for  each  level  of  Inspection. 


EVENT  A. LEVEL. INSPECTION 


DESCRIPTION  - This  event  represents  the  performance  of  an  A-level 
inspection.  The  constants  which  define  the  probability  of  detection 
are  passed  to  the  routine  EXAMINE  which  handles  the  actual  inspection. 

If  a defect  is  detected  on  one  of  the  ten  high  time  aircraft,  the  time  is 
stored  in  array  f!RDD,  This  event  can  be  scheduled  in  the  routine 
INSPECTION. SCHEDULER  and  the  events  A. LEVEL. INSPECTION,  B. LEVEL. INSPECTION, 
C. LEVEL.. INSPECTION,  and  D. LEVEL. INSPECTION. 

Local  Variables 


FOUND  - This  real  variable  Is  used  to  indicate  if  a defect  is  found 
and  the  type  of  defect. 

IDA  - This  integer  variable  is  the  indentification  number  of  the  aircraft 
being  inspected. 
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EVENT  B. LEVEL. INSPECTION 


DESCRIPTION  - This  event  represents  the  performance  of  a B-level  Inspection, 
the  Inspection  Is  not  done,  If  the  B-level  Inspection  Is  within  one-half  a 
B-level  interval  from  the  next  C-level  or  D-level  Inspection.  The  constants 
which  define  the  probability  of  detection  equation  at  the  B-level  are  passed 
to  the  routine  EXAMINE,  which  determines  If  the  defect  Is  found.  The 
Inspection  Interval  for  each  Inspection  level  Is  not  necessarily  an  even 
multiple  of  all  lower  level  intervals.  If  an  A-level  Inspection  is  scheduled, 
It  Is  cancelled  and  rescheduled  at  present  time,  plus  one  A-level  Interval 
later.  If  a defect  Is  detected  on  one  of  the  ten  high  time  aircraft,  the 
time  Is  stored  In  array  MRDD.  This  event  can  be  scheduled  In  the  routine 
INSPECTION. SCHEDULER  and  events  B. LEVEL. INSPECTION,  C. LEVEL. INSPECTION , 
and  D. LEVEL. INSPECTION. 

Local  Variables 

FOUND  - This  real  variable  is  used  to  indicate  if  a defect  Is  found  and 
the  type  of  defect. 

IDB  - This  Integer  variable  Is  the  identification  number  of  the  aircraft 
being  Inspected. 
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EVENT  C , LEVEL . INSPECTION 


DESCRIPTION  - This  event  represents  the  performance  of  a C-level 
Inspection.  The  Inspection  is  not  done,  if  the  C-level  Inspection  is 
within  one-half  a C-level  interval  of  the  next  D-level  inspection. 

If  there  is  either  an  A-level  or  a B-level  inspection  currently 
scheduled,  it  is  cancelled  and  rescheduled  at  one  A-level  or  B-level 
interval,  respectively,  later.  The  constants  which  define  the 
probability  of  defect  detection  equation  at  the  C-level  are  passed 
to  the  routine  EXAMINE,  which  determines  if  the  defect  is  found.  The 
C-level  inspection  uses  both  exploratory  and  non-expl oratory  probabil- 
ity of  detection  curves.  The  non-exploratory  curve  goes  into  effect 
after  one  crack  has  been  detected.  If  a defect  is  detected  on  one  of 
the  ten  high  time  aircraft,  the  time  is  stored  in  array  MRDD.  This 
event  can  be  scheduled  in  routine  INSPECTION. SCHEDULER  and  events 
C. LEVEL. INSPECTION  and  D. LEVEL. INSPECTION. 

Local  Variables 


FOUND  - This  real  variable  is  used  to  indicate  if  a defect  is  found 
and  the  type  of  defect. 

IDC  - This  integer  variable  is  the  identification  number  of  the  aircraft 
EeTng  inspected. 
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EVENT  D, LEVEL. INSPECTION 


DESCRIPTION  - This  event  represents  the  performance  of  a D-level  Inspec- 
tion!^ IT  there  are  any  lower  level  Inspections  currently  scheduled, 
they  are  cancelled  and  rescheduled  at  one  inspection  Interval  later.  A 
determination  Is  made  whether  an  aircraft  receives  an  external  or 
Internal  Inspection  dependent  on  the  current  value  of  the  variable 
SAMPLING.  The  constants  which  define  the  probability  of  defect  detection 
equation  at  the  D Internal  and  external  levels  are  passed  to  the  routine 
EXAMINE,  which  determines  if  a defect  is  found.  Explanatory  and  non- 
explanatory  probability  of  crack  detection  curves  are  used.  Non-explan- 
atory  crack  detection  curves  go  into  effect  after  a crack  has  been 
detected.  The  defect  histories  of  the  ten-high  time  aircraft  are  now 
examined.  If  all  of  the  ten-high  time  aircraft  time  goes  one  D-level 
interval  without  a defect  being  detected  at  any  inspection  level,  the  C- 
level  and  D-level  Intervals  are  increased  and  the  value  of  SAMPLING  is 
changed. 

The  changes  in  the  C and  D intervals  and  SAMPLING  are  dependent  on  the 
current  value  of  the  D-level  Interval.  If  the  current  D-level  Interval 
is  less  than  or  equal  to  12,000  hours,  the  following  relations  are  used: 


ABCD ( 3 ) 

NEW 

ABCD(3)  * 
OLD 

FREQ. DECREASE  /2  +1 

ABCD(4) 

NEW 

ABCD(4)  * 
OLD 

FREQ. DECREASE  *2  +1 

SAMPLING  = SAMPLING/(FREQ. DECREASE  *2  +1 )+l 
NEW  OLD 


If  the  current  D-level  interval  is  greater  than  12,000  hours,  the 
following  relations  are  used: 

ABCD(3)  = ABCD ( 3 ) * (FREQ. DECREASE  +1) 

NEW  OLD 

ABCD (4)  = ABCD(4)  * (FREQ. DECREASE  +1 ) 

NEW  OLD 

SAMPLING  = SAMPLING  +1 
NEW  OLD 

The  D-level  inspection  interval  Is  limited  to  32,000  hours  and  the  C- 
Level  interval  is  limited  to  32,000  times  the  ratio  of  the  initial  C- 
interval  to  the  initial  D-interval.  The  D-level  inspection  interval  for 
all  aircraft  is  rescheduled  at  the  new  value.  No  changes  are  permitted 
if  a modification  is  pending.  If  a service  modification  or  retrofit 
modification  is  pending,  It  is  done  at  the  D-level.  If  the  D-level 
interval  has  been  decreased  prior  to  a modification,  then  it  is  restored 
to  its  value  prior  to  the  decrease  after  the  modification. 
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A check  Is  done  In  the  routine  to  replace  the  ten  high  time  aircraft 
after  they  retire.  The  ten  high  time  aircraft  are  replaced  by  ten 
aircraft  whose  next  D-level  Inspection  occurs  closest  to  one  D-level 
interval  from  the  last  Increase.  After  this  point  the  last  ten  aircraft 
in  the  fleet  are  taken  as  the  high  time  aircraft. 

If  a defect  is  found  and  a modification  is  not  done,  a repair  is  scheduled. 
This  event  can  be  scheduled  in  the  events  ENTER. SERVICE,  D. LEVEL. 
INSPECTION,  T. INSPECTION. INCREASE  and  INCREASE. INSPECTION. FREQUENCY. 

Local  Variables 


DECC  - This  real  variable  is  the  factor  the  C-interval  is  multiplied  by 
for  an  interval  increase  when  old  D-interval  is  less  than  12,000  hours. 

DECD  - This  real  variable  is  the  factor  the  D-interval  is  multiplied  by 
for  an  interval  increase  when  old  D-interval  is  less  than  12,000  hours. 

FOUND  - This  real  variable  is  used  to  indicate  if  a defect  is  found  and 
the  type  of  defect. 

IDD  - This  integer  variable  is  the  identification  number  of  the  aircraft 
being  inspected. 

N - This  integer  variable  contains  the  identification  numbers  of  aircraft 
that  receive  internal  D-level  inspections. 

SMALL  - This  real  variable  is  used  in  determining  which  ten  aircraft 
replace  the  ten  high  time  aircraft  after  they  retire. 

MAJOR  LOGIC  STEPS 

1.  Schedule  the  next  D. LEVEL. INSPECTION. 

2.  Cancel  and  reschedule  lower  level  inspections  one  interval 
from  the  present  D-level  inspection. 

3.  Determine  if  internal  or  external  inspection. 

4.  Call  EXAMINE. 

5.  Check  high  time  aircraft  for  defect. 

6.  Increase  inspection  intervals  if  no  modification  pending 
and  no  defect  on  high  time  aircraft. 

7.  Reschedule  D. LEVEL. INSPECTION  if  step  6 is  true. 

8.  Determine  aircraft  that  will  replace  ten  high  time  aircraft 
after  they  retire. 

9.  If  modification  pending: 

a.  Call  INSTALL. MODIFICATION. 

b.  Reschedule  D. LEVEL. INSPECTION. 

c . Return . 

10.  If  defect  found,  schedule  REPAIR. 
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ROUTINE  EXAMINE 


DESCRIPTION  - This  routine  performs  the  numerical  comparison  which 
determines  whether  a defect  Is  detected  at  each  level  of  Inspection. 

This  routine  receives  constants  from  calling  routines  which  determine 
the  probability  of  defect  detection.  The  calculation  of  the  probability 
of  corrosion  detection  and  the  probability  of  crack  detection  that  is 
exponential  is  done  by  routine  PODD.  For  crack  lengths  that  are  not  on 
the  exponential  curve  the  probability  of  detection  Is  determined  In 
routine  EXAMINE  for  all  Inspection  levels.  Internal  defects  can  only  be 
detected  by  Internal  Inspection.  Cracks  that  originate  Internally  and 
become  external  can  be  detected  by  external  inspection,  but  use  only  the 
external  length  of  the  crack.  If  a defect  is  detected,  a repair  is 
scheduled  if  the  Inspection  level  is  A,  B,  or  C.  This  routine  can  be 
called  from  events  A. LEVEL. INSPECTION,  B. LEVEL. INSPECTION,  C. LEVEL. 
INSPECTION,  and  D. LEVEL. INSPECTION. 

Local  Variables 


AREA  - This  real  variable  is  the  calculated  area  of  the  corrosion 
defect. 

CCL1  - This  real  variable  Is  the  first  critical  crack  length. 

CCL2  - This  real  variable  is  the  second  critical  crack  length. 

CCL3  - This  real  variable  Is  the  third  critical  crack  length. 

a - This  real  variable  is  the  calculated  length  of  the  fatigue  crack. 

FOUND  - This  real  variable  is  set  equal  to  two  whenever  corrosion  is 
detected  and  to  one  whenever  a crack  is  detected. 

Ml , M2 , M3 , M4  - These  real  variables  are  the  individual  element  crack 
growth  rates. 

N - This  integer  variable  identifies  the  inspection  level  and  is  passed 
Fy  the  calling  event. 

PL  - This  real  variable  is  the  probability  of  defect  detection. 

TAC  - This  real  variable  is  the  simulation  time  of  corrosion  initiation. 

TA1  - This  real  variable  is  the  simulation  time  of  the  first  crack 
initiation. 

TA2-  This  real  variable  is  the  simulation  time  of  the  second  crack 
TnTtiation. 

XA,  XL,  VA,  ZA,  ZL  - These  real  variables  are  the  probability  of  defect 
detection  constants  for  the  exponential  curve  and  are  passed  by  the 
calling  routines. 
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1 - This  alpha  variable  Is  the  level  of  inspection  being  performed. 
MAJOR  LOGIC  STEPS 


1.  If  corrosion  exists,  and  inspection  Is  Internal,  or  corrosion 
is  external: 

a.  Calculate  corrosion  area. 

b.  Call  PODD  to  determine  probability  of  detection. 

c.  Draw  random  number  and  compare  to  probability  of 
detection.  If  random  number  less  than  or  equal  to 
probability  of  detection,  then  corrosion  is  found. 

2.  If  first  crack  exists,  and  inspection  is  internal  or  crack 
is  external : 

a.  Calculate  crack  length. 

b.  If  crack  originated  internally  and  inspection  is  external, 
subtract  internal  portion  of  crack  from  crack  length. 

c.  Determine  probability  of  detection  for  crack,  if  crack 
length  is  on  non-exponential  part  of  detection  curve. 
Probability  of  crack  detection  is  dependent  on 
inspection  level  and  crack  length. 

d.  If  crack  length  on  exponential  part  of  detection  curve, 
call  PODD  for  probability  of  detection. 

e.  Draw  random  number  and  compare  to  probability  of 
detection.  If  random  number  loss  than  or  equal  to 
probability  of  detection,  the  first  crack  is  found. 

f.  If  crack  found  and  crack  originated  internally  and 
external  inspection,  add  internal  portion  of  crack  back 
to  length  for  recording  of  total  crack  length. 

3.  If  second  crack  exists,  and  inspection  is  internal  or  crack 
is  external  (Same  as  steps  2A  - 2F). 

4.  If  defect  found,  schedule  REPAIR  if  Inspection  level  A,  13,  or 
C. 
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ROUTINE  POPP 


DESCRIPTION  - This  routine  computes  the  probability  of  detection  for 
corrosion  and  for  cracks  that  are  on  the  exponential  portion  of  the 
crack  detection  curve.  This  probability  Is  returned  to  the  calling 
routine.  PODD  Is  called  from  routine  EXAMINE  and  event  IMMEDIATE. 
FLEET. INSPECTION. 

Local  Variables 


J-  - This  real  variable  is  the  size  of  the  defect  under  consideration. 

RL  - This  real  variable  is  the  probability  of  detection  of  the  defect. 

X - This  real  variable  is  the  maximum  probability  of  detection  at  a given 
inspection  level . 

X - This  real  variable  is  an  empirically  determined  equation  constant 
for  each  level  of  inspection. 

1 - This  real  variable  is  the  minimum  defect  size  detectable  at  a given 
inspection  level. 
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ROUTINE  CANCEL . SCHEDULED. INSPECTIONS 


DESCRIPTION  - This  routine  cancels  all  scheduled  Inspections  below  the 
D- level  on  a given  aircraft.  Whenever  a defect  is  detected  and  repaired, 
it  is  assumed  that  all  other  defects  existing  on  that  particular  element 
are  also  repaired.  This  routine  is  called  at  this  time  to  cancel  all 
subsequent  inspections.  Also,  if  an  element  fails  or  an  aircraft  with 
existing  defects  is  retired,  this  routine  is  called  to  cancel  all 
scheduled  inspections.  This  routine  can  be  called  from  events  FAILURE, 
RETIRE. FROM. SERVICE,  and  REPAIR. 
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EVENT  REACH. FAIL. SAFE aST 


DESCRIPTION  - This  event  represents  the  time  when  the  residual  strength 
of  the  element  has  been  reduced  to  the  fill -safe  strength.  The  time  and 
aircraft  identification  numbers  are  saved  as  output.  The  calculations 
of  the  strength  reduction  Is  based  on  the  sum  of  the  two  crack  lengths, 
if  the  second  crack  Is  present  In  the  element.  This  event  can  be 
scheduled  in  events  1 .STRENGTH. REDUCTION  and  2. STRENGTH. REDUCT I ON. 

Local  Variables 

IDRFS  - This  Integer  variable  Is  the  Identification  number  of  the  aircraft 
being  processed. 
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EVENT  FAILURE 


DESCRIPTION  - This  event  represents  structural  failure.  When  this  event 
occurs,  the  aircraft  is  removed  from  the  active  fleet.  If  this  aircraft 
was  one  of  the  ten  high-time  aircraft  being  monitored  for  the  purpose  of 
increasing  inspection  intervals,  it  is  replaced  by  the  next  high-time 
aircraft  in  active  service.  Any  remaining  scheduled  events  are  cancelled 
and  their  event  routines  destroyed.  The  events  INCREASE. INSPECTION. 
FREQUENCY  and  routine  COMP. RISK  is  called.  This  event  can  be  scheduled 
in  events  1 .STRENGTH. REDUCTION  and  2.  STRENGTH. REDUCTION. 

Local  Variables 

HOLD  This  integer  variable  serves  as  an  intermediate  storage  for  aircraft 
identification  when  replacing  one  of  the  ten  high-time  aircrafts. 

IDFA  - This  integer  variable  is  the  identification  number  of  the  aircraft 
under  consideration. 

S - This  real  variable  is  the  residual  strength  at  failure. 


i 
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EVENT  RETIRE. FROM. SERVICE 


DESCRIPTION  - This  event  represents  the  retirement  of  an  aircraft  from 
active  service.  If  a crack  exists  the  routine  COMP. RISK  is  called.  The 
aircraft  being  retired  is  removed  from  the  set  ACTIVE. FLEET  and  filed  In 
the  set  FLEET. RETIRED.  All  remaining  scheduled  events  for  this  aircraft 
are  cancelled  and  the  event  notices  destroyed.  This  event  Is  scheduled 
in  the  event  ENTER. SERVICE. 

Local  Variables 


IDRET  - This  integer  variable  Is  the  Identification  number  of  the 
aircraft  being  retired  from  service. 


- 81  - 


EVENT  REPAIR 


DESCRIPTION  - This  event  represents  the  structural  repair  of  an  element. 
Tf  the" events  FAILURE  and  REACH. FAIL. SAFE. LQT  are  scheduled,  they  are 
cancelled.  It  Is  assumed  that  all  existing  defects  are  repaired  and 
that  new  times  to  defect  occurrences  are  determined  In  the  same  manner 
as  when  the  aircraft  entered  service.  The  size  and  number  of  all 
existing  cracks  are  compared  with  the  Inspection  interval  decrease 
criteria. 

The  inspection  interval  decrease  criteria  differentiates  between  large 
and  small  cracks.  Large  crack  criteria  is  a crack  greater  than  the 
fail-safe  length  or  if  the  sum  of  the  crack  lengths  plus  crack  growth 
rate  times  inspection  interval  Is  greater  than  (SU-1 .2*(DLL+1 )-l )*F$AF. 
LGT/ ( SU-SF ) . The  inspection  interval  is  the  current  D-level  interval 
if  the  crack  is  found  in  an  external  inspection.  Otherwise,  the 
inspection  interval  is  the  current  D-level  interval  times  the  current 
SAMPLING.  If  a large  crack  is  found,  an  IMMEDIATE. FLEET. INSPECTION 
and  INCREASE. INSPECTION. FREQUENCY  is  scheduled. 

If  the  special  fleet  inspection  detects  a crack,  the  inspection  Interval 
can  be  reduced  again  if  the  crack  detected  during  the  special 
inspection  meets  the  criteria  just  discussed.  Small  cracks  are  those 
not  determined  by  the  criteria  to  be  large.  If  the  sum  of  the  cracks 
found  in  the  entire  fleet  times  the  fail-safe  crack  length  is  greater 
than  one-fifth  of  the  fleet  size,  the  inspection  interval  is  reduced. 

If  a modification  is  pending  on  the  aircraft  being  repaired,  the 
modification  is  installed  at  this  time.  A decision  on  modification,  due 
to  the  current  defect,  is  scheduled.  The  routine  COMP. RISK  is  called. 
This  event  can  be  scheduled  in  routine  EXAMINE  and  events  D. LEVEL. 
INSPECTION  and  IMMEDIATE. FLEET. INSPECTION. 

Local  Variables 


AAFL  - This  real  variable  is  the  actual  average  fatigue  life  of  the  element. 
CCL1 , CCL2,  CCL3  - These  real  variables  are  the  critical  crack  lengths. 

CL  - This  real  variable  is  the  crack  length. 

FIRST. LIFE  - This  real  variable  is  the  time  in  flight  hours  to  first 
crack  initiation. 

HOURS. TO. CORROSION  - This  real  variable  is  the  time  in  flight  hours  to 
corros ion  initTatTon . 
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I PREP  - This  integer  variable  Is  the  aircraft  Identification  number. 


MAX.CRK  - This  real  variable  Is  the  sum  of  the  crack  lengths  In  the 
element. 

POT.CRK  - This  real  variable  is  the  sum  of  the  maximum  crack  length  and 
the  length  that  the  crack  will  gain  In  the  next  D-level  Inspection 
interval. 

RST  - This  real  variable  is  the  remaining  service  time  of  the  aircraft. 
SDL  - This  real  variable  is  1.2*  DLL  - 1. 

SECOND. LIFE  - This  real  variable  Is  the  time  In  flight  hours  to  the 
second  crack  Initiation. 

STR.P.ED  - This  real  variable  is  the  element  strength  reduction  because 
of  all  existing  cracks. 

TAC  - This  real  variable  Is  simulation  time  to  corrosion  initiation. 

TA1 , TA2  - These  real  variables  are  simulation  time  to  the  first  and 
second  crack  initiation. 

MAJOR  LOGIC  STEPS 

1.  Cancel  REACH. FAIL. SAFE. LGT  and  FAILURE,  if  scheduled. 

2.  Cancel  CORROSION  if  corrosion  exists. 

3.  Call  PREDICT. CORROSION. 

4.  Schedule  CORROSION  if  It  occurs  within  remaining 
service  life. 

5.  Calculate  first  and  second  crack  lengths. 

6.  Calculate  fleet  strength  reduction  and  potential 
crdck  growth. 

7.  If  crack  large,  schedule  an  IMMEDIATE. FLEET. INSPECTION  and 
an  INCREASE. INSPECTION. FREQUENCY. 

8.  If  crack  small  and  fleet  strength  reduction  criteria  met, 
schedule  an  INCREASE. INSPECTION. FREQUENCY. 

9.  Install  modification  if  pending,  call  INSTALL. MODIFICATION. 

10.  For  fatigue  test  failure  modification  cancel  T. INSPECTION. 
INCREASE  if  scheduled. 

11.  If  crack  exists,  call  COMP. RISK. 

12.  Schedule  DECISION. ON. MOD. 

13.  Cancel  and  reschedule  fatigue  crack  initiation. 

14.  If  inspection  scheduled,  call  CANCEL. SCHEDULED. INSPECTIONS. 
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EVENT  T. INSPECTION. INCREASE 


DESCRIPTION  - This  represents  an  Inspection  interval  decrease  for  a 
particular'  aircraft  pending  a structural  modification  because  of  a 
fatigue  test  failure.  The  D-Interval  is  reduced  by  the  Input  quantity 
T.FREQ.CHG.  The  D-level  inspection  is  rescheduled  for  aircraft  with 
modification  pending.  This  event  is  scheduled  in  event  T. IMPLEMENT. MOD. 

Local  Variables 

I DTI  - This  integer  variable  is  the  identification  number  of  the  aircraft 
under  consideration. 


EVENT  INCREASE. INSPECTION. FREQUENCY 


DESCRIPTION  - This  event  represents  a fleet  wide  decrease  In  the  D-level 
TnspectTon”1nterval  and  a change  In  SAMPLING.  The  D-level  Inspection 
Interval  Is  reduced  by  the  Input  quantity  S. FREQ. CHG. SAMPLING  Is  set 
equal  to  one,  If  the  crack  is  detected  In  a D-level  Internal  Inspection. 

For  SAMPLING  equal  to  one,  every  aircraft  Is  Inspected  Internally  at  the 
D-level.  If  the  crack  Is  not  found  in  an  Internal  D-leve"'  Inspection 
SAMPLING  is  changed  by  dividing  Its  current  value  by  S. FREQ. CHG.  The  D- 
level  Inspection  Interval  for  the  entire  fleet  Is  then  rescheduled. 

This  event  can  be  scheduled  in  events  REPAIR  and  IMMEDIATE. FLEET. INSPECTION. 
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EVENT  IMMEDIATE. FLEET. INSPECTION 


DESCRIPTION  - This  event  represents  an  Immeoiate  fleet-wide  inspection 
caused  by  finding  a large  crack.  Existing  crack  lengths  and  corrosion 
areas  are  calculated  along  with  the  associated  probabilities  of  detec- 
tion. As  in  the  scheduled  inspections,  these  probabilities  are  compared 
with  a random  number  to  determine  whether  or  not  the  defect  is  detected. 
If  any  cracks  are  found  and  the  inspection  Intervals  had  not  previously 
been  reduced  an  INCREASE. INSPECTION. FREQUENCY  is  scheduled.  If  any 
defects  are  found  a REPAIR  is  scheduled.  These  defects  will  not  cause 
any  additional  fleet-wide  inspections.  This  event  is  scheduled  in  event 
REPAIR. 

Local  Variables 


AREA  - This  real  variable  is  the  area  of  an  existing  corrosion  defect. 

CCL1 , CCL2,  CCL3  - These  real  variables  are  the  critical  crack  lengths. 

CK.FIND  - This  real  variable  is  set  equal  to  one  if  a crack  is  found  in 
the  special  inspection. 

CL  - This  real  variable  is  the  length  of  an  existing  crack. 

FOUND  - This  real  variable  indicates  if  a defect  is  found  and  the  type 
of  defect. 

Ml , M2 , M3 , M4  - These  real  variables  are  the  four  crack  growth  rates. 

PL  - This  real  variable  is  the  probability  of  detection  of  an  existing 

defect. 

TAC  - This  real  variable  is  the  time  of  Initiation  of  an  existing 
corrosion  defect. 

TA1 , TA2  - These  real  variables  are  the  times  to  initiation  of  the  first 
and  second  cracks. 

MAJOR  LOGIC  STEPS 

1 . If  corrosion  exists: 

a.  Calculate  area. 

b.  Determine  probability  of  detection,  call  PODD. 

c.  Compare  with  random  number  to  determine  if  detected. 

2.  If  first  crack  exists: 

a.  Calculate  length. 

b.  Determine  probability  of  detection,  call  PODD  if 
on  exponential  part  of  detection  curve. 

c.  Compare  with  random  number  to  determine  if  detected. 

3.  If  second  crack  exists: 

(Same  as  2A  - 2C) . 
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4.  If  high  time  aircraft  have  defects,  record  time. 

5.  If  defects  found,  schedule  REPAIR. 

6.  If  crack  found  and  Inspection  Interval  not  previously 
reduced,  schedule  an  INCREASE. INSPECTION. FREQUENCY. 
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EVENT  DEC I SION. ON. MOD 


This  event  makes  the  decision  on  whether  or  not  to  develop  a structural 
modification  because  of  service  experience.  The  decision  to  develop  a 
modification  is  made  by  comparing  the  cost  per  flight  hours  of  the 
modification  with  the  repair  cost  per  flight  hours  plus  the  increased 
inspection  cost  per  flight  hours.  The  modification  cost  per  flight 
hours  is  found  by  dividing  the  total  fleet  modification  cost  by  the 
remaining  service  life  of  the  fleet.  The  repair  cost  per  flight  hours 
is  found  by  dividing  the  total  fleet  repair  costs  since  the  last  modi- 
fication by  the  fleet  flight  time  since  the  last  modification.  The 
increased  inspection  cost  per  flight  hour  is  found  by  dividing  the 
projected  increased  inspection  costs  by  the  remaining  service  life  of 
the  fleet.  SAMPLING  Is  included  in  the  Increased  inspection  costs 
calculation  since  cracks  found  in  Internal  D-level  inspections  set 
SAMPLING  equal  to  one,  while  SAMPLING  is  changed  by  S.FREQ.CHG  if  the 
crack  is  not  found  at  the  internal  D-level  inspection.  If  a decision  is 
made  to  modify,  the  event  IMPLEMENT. MODIFICATIuN  is  scheduled.  This 
event  is  scheduled  in  the  event  REPAIR. 

LOCAL  VARIABLES 


ACCUMULATED. HRS  - This  real  variable  is  the  total  fleet  time  since  the 
last  modification. 

ICPH  - This  real  variable  is  the  increased  inspection  cost  per  flight 
hour. 

MCPH  - This  real  variable  is  the  modification  cost  per  flight  hour. 

MD. COST  - This  real  variable  is  the  cost  of  installing  a modification  on 
a single  aircraft. 

MRFH  - This  real  variable  is  the  service  time  remaining  on  a particular 
aircraft  after  its  modification. 

NFTS  - This  real  variable  is  the  total  production  time  of  all  aircraft 
entering  service  after  the  second  production  rate  goes  into  effect. 

NPDL  - This  real  variable  is  the  number  of  aircraft  which  have  entered 
service. 

POST. MOD. HRS  - This  real  variable  is  the  total  fleet  service  time 
remain! ng'aTter  the  modification. 

RCPH  - This  real  variable  is  the  repair  cost  per  flight  hour. 

TOOLING  - This  real  variable  is  the  tooling  cost  in  the  development  of  a 
modification. 
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EVENT  IMPLEMENT . MODIFICATION 


DESCRIPTION  - i his  event  represents  the  development  of  a modification 
because  of  service  experience.  If  the  modification  Is  to  be  fatigue 
tested,  the  actual  average  fatigue  life  is  set  equal  to  the  original 
predicted  life  of  the  element  design.  Otherwise,  the  actual  average 
fatigue  life  Is  determined  by  calling  routine  REAL. LIFE.  The  modi- 
fication will  be  Installed  on  each  aircraft  at  the  next  D-level  Inspec- 
tion or  defect  repair.  The  D-level  inspection  interval  and  D-level 
internal  sampling  are  set  to  their  values  prior  to  the  interval  reduc- 
tion. This  event  can  be  scheduled  in  event  DECISION. ON. MOD. 

Local  Variables 


NS I G - This  real  variable  is  the  standard  deviation  of  the  ratio  distribu 
t ion  passed  to  routine  REAL. LIFE. 

NMU  - This  real  variable  is  the  mean  of  the  ratio  distribution  passed  to 
routine  REAL. LIFE. 


ROUTINE  DISPLAY. OUTPUT 


DESCRIPTION  - This  routine  prints  the  output  for  each  element.  This 
output  is  suppressed  If  the  long  list  option  Is  in  effect.  This  routine 
is  called  from  the  MAIN  program. 

Local  Variabl es 


All  the  local  variables  in  this  routine  are  used  to  temporarily  store 
output  values. 
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ROUTINE  SUMMARY 


DESCRIPTION  - This  routine  prints  the  summary  output  for  each  element 
type.  This  output  Is  suppressed  If  the  long  list  option  Is  In  effect. 
This  routine  Is  called  from  the  MAIN  program. 


Local  Variables 

All  the  local  variables  in  this  routine  are  used  to  temporarily  store 
output  values. 


ROUTINE  SUMMARY 


DESCRIPTION  - This  routine  prints  the  summary  output  for  each  element 
type.  This  output  is  suppressed  if  the  long  list  option  is  In  effect. 
This  routine  is  called  from  the  MAIN  program. 

Local  Variables 


All  the  local  variables  in  this  routine  are  used  to  temporarily  store 
output  values. 
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